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Research  Objectives 

The  broad  objectives  of  this  program  are  to  further  our  understanding  of  the 
intriiEie  and  extrinsic  properties  of  highly  transparent  solids  in  the  infrared 
region  of  the  spectrum.  The  identification  of  chemical  and  structural  imperfec¬ 
tions  which  occur  in  highly  transparent  insulators  and  semiconductors  are  in¬ 
vestigated  by  the  following  experimental  techniques:  infrared  optical  deriva¬ 
tive  and  light  scattering  techniques.  The  experimental  results  are  used  to 
develop  theoretical  models  for  the  intrinsic  and  extrinsic  interactions. 

Status  Report  Summary 

The  previous  informal  technical  report  covered  the  period  1  January  1980 
to  31  August  1980.  Since  this  interim  report  covers  the  period  1  January  1980 
to  31  March  1981,  seme  aspects  of  the  former  report  are  included  for  the  pur¬ 
poses  of  continuity. 

Our  unique  infrared  wavelength  modulated  spectrometer  is  capable  of  detect¬ 
ing  a  change  in  absorption  or  reflection  of  a  part  in  10  on  a  relatively  smooth 
background  in  the  spectral  region  from  0.2-2-20  Urns.  With  this  instrumentation 
we  have  for  the  first  time  measured  the  continuous  spectral  distribution  of  the 
absorption  in  highly  transparent  solids  at  levels  of  10  ^cm  A  detailed 
study  has  been  completed  on:  KC?.,  KBr,  CaF2>  LiF,  NaCi,  LaF^.  BaF^,  MgF0,  SrF0 
and  MgO  in  the  spectral  region  2.5-12.0  Jm.  Previously,  only  measurements  at 
discrete  wavelengths  by  thermocouple  laser  calorimetry  and  photo-acoustic  calori¬ 
metry  were  available.  The  richness  of  the  absorption  structures  observed  in  all 
these  crystals  enabled  us  to  identify  volume  and  surface  contaminants  and  allowed 
the  identification  of  the  type  of  bonding  of  the  particular  constituents  present. 
The  sensitivity  of  the  infrared  modulation  systems  indicates  that  it  is  capable 
of  detecting  a  fraction  of  a  monolayer  surface  species. 

Stress  in  chemical  vapor  deposition  layers  of  Si  on  sapphire  were  examined 
using  our  derivative  spectrometer  by  measuring  the  shift  and  polarization  depend¬ 
ence  of  the  3.4  eV  critical  point. 

Strain  in  ion  implanted  GaAs  was  readily  detected  by  our  techniques  bv  ob¬ 
serving  the  shift  of  the  critical  point. 
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Free-carr:ler  screening  of  hyperbolic  excicons  in  GaAs  was  observed. 

Deep  levels  in  Si  and  GaAs  were  observed  by  deep  level  derivative  absorp- 
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tion  spectroscopy  (DLDAS)  with  a  sensitivity  of  10  “-10  impurities/cm  pre¬ 
viously  only  obtained  by  electrical  transient  capacitance  techniques  (DLTS) . 

Oxides  on  the  compound  semiconductors  GaAs  and  HgCdTe  were  observed  in 

_  O  _ 

700  A  layers  using  Raman  scattering  and  optical  derivative  techniques;  OH” can 
be  observed  in  the  former  and  TeO^  in  the  latter  meterials.  In  addition,  SiH 

O 

was  observed  in  8A  layers  of  Si. 

An  initial  study  of  Rayleigh-Brillcuir.  scattering  from  single  and  poly¬ 
crystalline  alkali  halides  was  completed  indicating  that  dislocations,  grain 
boundaries  and  the  impurity-pinning  at  these  sites  are  possible  candidates  for 
increased  scattering  in  polycrystalline  samples 

A  series  of  infrared  transmitting  low  loss  fluoride  glasses  were  prepared 

using  reactive  atmosphere  processing  techniques.  In  addiLicn,  Rayleigh-Brillouin 

and  Raman  scattering  were  performed  on  these  glasses  indicated  that  the  HfF,- 

BaF0-ThF  system  scatters  more  than  the  ZrF  -ThF  -BaF„  svstem. 

1  4  4  4  2 

The  study  of  surface  polaritons  cn  n+  Silicon  and  Copper  was  completed 
showing  chat  such  studies  should  be  conducted  in  the  complex  plane  when  either 
the  damping  is  large  or  there  exists  iriterband  transitions  in  the  spectral  re¬ 
gion  studied. 

A  multipass  tandem  Fabry-Perot  interferometer  was  assembled  to  study  sur¬ 
face  excitations  on  semiconductors  as  well  as  the  largely  unexplored  region  of 
interface  surface  phonons  between  metal-conductor,  metal-oxide-serciconductor 
and  heterostructures  using  Rayleigh-Brillcuir.  lack-scattering. 

Infrared  Wavelength  Modulation  Spectroscopy  of  Insulators 

The  infrared  wavelength  modulation  specr -ometer  system  that  we  have  de¬ 
veloped  has  enabled  us  to  perform  unique  measurements  on  highly  transparent 
solids  at  levels  of  absorption  of  10  ^cm.  These  techniques  were  previously  em¬ 
ployed  on  KBr  and  KC2.  and  have  now  been  extended  to  other  crystalline  systems 
which  are  of  interest  in  light  guiding  applications.  A  detailed  study  has  been 
completed  for  CaF^,  LiF,  NaCi,  NaF,  LaF^,  3aF^,  MgF,,,  SrF^  and  MgO.  We  have 
measured  the  continuous  spectral  distribution  of  the  absorption  in  these  mate¬ 
rials  in  the  spectral  region  from  2:5  to  12.0  ym.  Previously,  only  measurements 


at  discrete  wavelengths  by  thermocouple  laser  calorimetry  and  photo-acoustic 
calorimetry  were  available.  Rich  and  varied  adsorption  structures  were  observed 
in  all  these  crystals  which  enable  us  to  identify  volume  and  surface  absorption. 

It  is  interesting  to  note  that  despite  the  diverse  chemistry  of  these  substances 
and  the  varied  crystal  growth  techniques  employed  in  the  preparation  of  these 
materials,  similar  dominant  bands  are  observed  in  most  samples  which  are  due  to 
physisorption,  volume  and  surface  chemisortpion.  However,  the  various  fine 
structures  observed  in  different  crystals  are  indicative  of  the  individual  char¬ 
acteristics  of  the  chemistry  of  the  crystal  preparation.  These  measurements 
which  were  performed  in  the  laboratory  and  dry  ambient  readily  show  the  physi¬ 
sorption  and  desorption  of  surface  contaminants  and  allow  the  identification  of 

the  type  of  bonding  of  the  particular  constituents  present.  The  surface  charac- 
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ter  of  some  of  the  absorption  at  levels  of  10  to  10  cm  indicates  the  im¬ 
portance  of  performing  measurements  in  the  same  ambient  when  comparing  low  level 
absorption  measurements  by  various  sensitive  techniques. 

The  sensitivity  of  the  infrared  modulation  techniques  indicates  that  a 
fraction  of  a  monolayer  of  surface  adsorbed  species  can  be  detected.  Experi¬ 
ments  in  controlled  ambients  on  well  characterized  substrates  using  this  tech¬ 
nique  can  be  rewarding  for  the  physics  and  chemistry  of  surfaces. 

This  work  has  been  presented  at  the  Laser  Damage  Symposium,  1980,  at  Boulder, 
Colorado.  A  paper  will  appear  in  the  NBS  Proceedings  of  this  Conference.  A 
copy  of  the  submitted  manuscript  is  enclosed  in  the  Appendix  of  this  report. 

The  successful  use  of  our  infrared  optical  derivative  system  in  the  study 
of  bulk  and  surface  impurities  in  insulators  at  levels  of  10  "^cra  indicates  that 
much  remains  to  be  done  in  the  area  of  laser  windows.  It  is  clear  that  our 
-modulation  measurements  should  be  used  in  conjunction  with  a  systematic  study 
with  crystal  synthesis  groups. 

In  the  latter  part  of  this  report  period  our  attention  was  directed  toward 
the  use  of  our  infrared  A-modulation  techniques  as  well  as  Ravleigh-Brillouin 
scattering  to  semiconductor  problems.  We  have  demonstrated  that  A-modulation 
can  be  used  to  study  deep  levels  in  GaAs  And  silicon  at  levels  of  sensitivitv 
formerly  only  obtainable  using  DLTS  electric  measurements.  In  addition,  we 
have  shown  that  thin  surface  oxides  can  be  studied  by  A-modulation  and  optical 
scattering  techniques. 
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Wavelength  Spectroscopy  of  Semiconductors 

In  the  growth  of  semiconductors  by  various  expitaxial  techniques  and  the 
doping  by  icn  implantation,  an  important  technological  problem  is  a  means  of 
detecting  stress  in  the  growth  layers.  We  have  used  our  wavelength  modulation 
techniques  to  determine  the  shift  of  various  band  structure  critical  points  as 
a  function  of  doping.  It  has  been  possible  to  observe  the  stress  in  chemical 
vapor  deposition  layers  (CVD)  of  Si  on  sapphire  by  examining  the  shift  of  the 
3.4  eV  critical  point.  In  addition  a  polarization  dependence  of  the  shift  was 

O 

measured  indicating  an  anisotropic  strain  exists  in  4700  to  1900  A  layers  on  Si. 

Strain  in  ion  implanted  GaAs  was  readily  observed  by  our  techniques.  The 

shift  in  the  L^-L^  critical  point  due  to  various  ion  implanted  species  such  as 

Be,  Sb,  S,  In  and  double  implanted  Sb  +  Be  were  observed;  the  implanting  flux 

13  2 

was  of  the  order  of  10  /cm  .  It  is  surprising  that  significant  structural 
changes  were  observed  for  such  low  levels  of  implanting  and  attests  to  the  sen¬ 
sitivity  of  cur  techniques.  Compared  to  unimplanted  GaAs,  positive  and  nega¬ 
tive  shifts  of  the  energy  of  the  critical  point  indicates  that  we  are  able  to 
distinguish  contraction  and  expansion  of  the  lattice.  Previous  attempts  to 
study  strain  effects  in  semiconductors  looked  for  changes  in  carrier  mobility 
in  implanted  material  proved  futile.  Since  transport  measurements  involve  aver¬ 
ages  over  the  distribution  of  Bloch  functions,  the  local  character  of  strain  is 
not  immediately  obvious.  Since  our  wavelength  modulation  reflectivity  measure¬ 
ments  samples  the  material  within  a  penetration  depth  for  the  radiation,  in  the 

O 

present  case  **  5CA,  one  can  "view"  local  effects  of  strain  on  the  band  structure 
of  the  semiconductor. 

We  have  begun  studies  of  the  E^  and  E^  +  critical  point  in  GaAs  as  a 
function  of  doping.  Since  this  van  Hove  singularity  has  a  hyperbolic  exciton 
associated  with  it,  it  is  of  interest  to  study  the  effects  of  free-charge  screen¬ 
ing  on  the  final  state  interaction.  This  study  has  shown  that  the  reflectivity 
as  a  function  of  increasing  carrier  concentration,  shows  little  broadening,  but 
a  definite  decrease  in  magnitude.  These  effects  are  indicative  of  screening 
effects  on  the  hyperbolic  exciton.  Although  there  have  been  many  experimental 
and  theoretical  studies  to  search  for  free-carrier  screening  on  .direct  and  in¬ 
direct  excitons,  to  our  knowledge  none  have  concerned  the  hyperbolic  exciton  of 
this  critical  point.  This  is  particularly  surprising  in  view  of  the  extensive 
studies  of  band  structure  of  semiconductors  by  electro,  piezo,  thermo-modulation 


techniques.  It  is  planned  to  continue  these  studies  foi  n  and  p-tvpe  doped 
materials.  In  addition,  we  plan  to  develop  the  theory  of  screening  for  this  ex- 
citon.  The  study  of  screening  effects  for  such  final  state  interactions  is  par- 
ticulary  attractive  since  one  can  avoid  the  ancillarv  effects  of  screening  such 
as  bam  "ailing  or  Burstein  shifts  which  can  mask  screening  effects  for  direct 
and  indirect  excitons.  Knowledge  of  screening  effects  will  become  increasingly 
important  technologically  as  very  large  scale  integration  of  heavily  doped  de¬ 
vices  develops. 


Deep  Level  Derivative  Absorption  Spectroscopy  (DLDAS)  on  GaAs  and  Si 

Our  understanding  of  point  defects,  such  as  vacancies  and  chemical  impuri¬ 
ties  in  semiconductors  which  give  rise  to  substantial  rearrangements  of  elec¬ 
tronic  density  and  atomic  positions,  are  not  well  understood  since  they  are 
not  well  described  by  effective  mass  or  linear  response  theory.  In  recent 
years  the  emergence  of  transient  capacitance  techniques  have  stimulated  quan¬ 
titative  studies  of  deep  levels  and  motivated  theoretical  studies.  Despite  the 
immense  variation  in  these  junction  techniques  which  differ  in  technical  de¬ 
tails,  the  end  results  one  desires  from  these  measurements,  is  the  optical  emis¬ 
sion  or  absorption  crossections  for  either  electrons  or  holes  which  are  pro¬ 
portional  to  the  absorption  coefficient.  In  general,  excited  states  are  not 
observed  in  these  techniques  and  thresholds  for  the  transitio-  s  to  bands  are  not 
easily  determined.  Direct  absorption  measurements  could  yield  the  quantities 
of  interest,  but  at  levels  of  sensitivity  of  DLTS  techniques,  on  the  order  of 
10^-10^  /cm\  they  can  not  be  employed.  However  our  optical  derivative  tech¬ 
niques  have  the  requisite  sensitivity.  In  addition  it  is  not  necessary  to  make 
electrical  contacts  to  the  samples  nor  are  we  restricted  to  specific  sample  re- 
sisitivies . 

We  have  performed  wavelength  derivative  measurements  on  high  purity  triple 

zone  refined  Si  and  have  observed  numerous  absorption  bands  between  the  band 
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edge  and  the  multiphonon  tail  at  levels  of  10  cm  .  The  identification  of 
these  levels  have  not  been  definitively  made;  some  of  them  correspond  to  Si 
vacancies  and  interstitial  energies.  We  plan  further  measurements  on  Si  as  a 
function  of  temperature.  Aside  from  the  above  bulk  leveLs  we  have  also  ob¬ 
served  an  ahsorption  band  at  4.4  um  which  is  present  in  Si  whose  strength  is 
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independent  of  sample  thickness.  This  band  seems  to  be  due  to  an  Si-H  vibra¬ 
tional  band;  it  is  particularly  strong  when  Si  is  etched  with  HF,  however  it 
is  present  regardless  of  the  surface  preparation  employed. 

The  properties  of  high  resistivity  GaAs  continues  to  be  of  technological 
importance  since  it  is  a  substrate  or  buffer  layer  in  molecular  beam  and  other 
vaporous  growth  techniques.  The  aim  of  electrical  characterization  of  insulating 
materials  is  to  determine  the  concentration  of  shallow  donor  or  acceptors  to  be 
compensated  by  deep  levels.  To  this  day  it  is  not  clear  whether  Cr  or  0  levels 
used  for  compensation  give  rise  to  multiple  levels  or  exist  as  clusters. 

We  have  initiated  studies  by  optical  derivative  techniquest  to  study  semi- 
insulating  GaAs.  In  CrrGaAs  samples  20  mils  thick,  we  have  easily  resolved  a 
number  of  absorption  bands  between  0.8  and  1.1  eV.  In  addition,  absorption  bands 

were  observed  in  undoped  crystals  with  absorption  coefficients  at  levels  of 
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10  cm  .  Most  previous  studies  of  such  material  have  only  been  studied  by  OTSC 
techniques  (optical  thermally  stimulated  currents).  We  have  initiated  such  elec¬ 
trical  studies  in  conjunction  with  our  optical  derivative  techniques.  Simultaneous 
studies  using  optical  and  thermal  stimulation  can  be  profitable  since  each  yield 
slightly  different  physical  parameters  for  the  deep  levels. 

In  addition,  we  have  initiated  high  resolution  Fourier-transforn  photocon¬ 
ductivity  measurements  to  look  for  the  sample  differences  in  the  hvdrogen-like 
spectra  of  shallow  donor  or  acceptor  species  in  MBS  prepared  layers. 

The  tripartite  approach  or  the  same  sample  will  give  a  global  assessment 
of  semi-insulating  GaAs.  The  purpose  of  such  an  assessment  is  to  determine  the 
residual  concentration  of  shallow  levels  15^-N^  ar>d  t*"'e  concentration  °f  compen¬ 
sating  deep  levels. 


Raman  Scattering  and  Optical  Derivative  Measurements  of 
Oxides  on  HgCdTe,  GaAs  and  Si 

The  oxides  that  are  formed  on  compound  semiconductors  are  of  fundamental 
and  technological  interest.  The  sites  of  chemisorption  on  these  materials  have 
been  a  matter  of  discussion  for  some  time.  The  successful  fabrication  of  MOS 
type  devices  as  well  as  the  stabilization  of  surface  states  depends  upon  the 
preparation  of  appropriate  oxides.  The  chemical,  phvsical  and  electrical  pro- 
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pertius  of  oxides  on  C-aAs  and  HgCdTe  are  noc  sufficiently  understood  or  con¬ 
trolled.  The  use  of  microprobe ,  Auger  and  other  commonly  used  surface  analysis 
methods  yield  only  atomic  abundance.  However,  Raman  spectroscopy  can  provide 
direct  information  on  the  composition  and  structure  of  surface  oxides. 

O 

We  have  performed  Raman  measurements  on  a  sample  of  700  A  oxide  grown  on 
HgCdTe  obtained  from  the  Santa  Barbara  Research  Center.  In  this  work  a  number 
of  bands  were  observed  in  the  region  100  to  400  cm  To  identify  these  oxides, 
a  series  of  Raman  experiments  were  performed  on  crystalline  CdO  and  Te0o  since 
there  is  a  paucity  of  literature  of  the  modes  of  these  oxides  which  are  suspected 
to  occur  on  HgCdTe  surfaces.  This  work  has  allowed  us  to  deconvoiute  some  or  the 
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structures  observed  in  the  700  A  layer.  The  dominant  bands  observed  in  these 
oxides  seem  to  be  due  to  oxygen  bound  to  Te-sites.  We  have  also  been  able  to 
observe  infrared  active  bands  in  these  oxides  using  infrared  derivative  tech¬ 
niques. 

In  examining  the  oxides  grown  on  GaAs  substrates,  it  has  been  possible  to 
detect  OH  modes  in  0.5  ym  layers.  We  have  previously  shown,  using  our  infrared 
modulation  system  to  study  the  9.0  ym  SiO  vibrational  band  in  °..l  1  avers  "nas¬ 
cent  oxide  on  Si  form  in  normal  laboratorv  -ambient.  In  these  experiments,  5i 
substrates  were  etched  with  HF  and  the  growth  of  the  oxide  was  monitored  by  ob¬ 
serving  the  growth  of  the  9,0  Um  band.  As  we  have  previously  indicated  in  this 
report5we  also  observed  strong  SiK  bands  when  HF  etch  was  used  as  well  as  when 
other  surface  preparations  were  employed. 

It  is  clear  that  since  our  methods  can  detect  a  fraction  of  a  monolayer  of 
oxide,  if  we  were  to  perform  such  experiments  in  UHV  we  could  answer  the  vexing 
question  involved  in  the  passevation  of  GaAs  surfaces  as  to  what  is  the  initial 
site  of  oxidation  on  GaAs.  By  observing  whether  tiioO-  orAs-iOg  hands  are  initial!', 
formed  in  UHV  oxidation,  an  important  contribution  to  this  field  can  be  made.  Wo 
have  given  some  consideration  to  performing  such  UHV  optical  experiments.  The 
major  components  of  UHV  pumping  and  dewar  system  have  been  assembled.  However 
some  auxilltary  equipment  such  as  a  quadrupoLe  mass-spectromemter  is  not  avail¬ 
able  to  complete  such  a  system. 
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Ray leigh- Bril loui n  Scattering 

We  have  assembled  a  high  contrast  ratio  live-pass  Fabrv-Perot  interfero¬ 
meter.  In  addition,  a  tandem  Brillouin  spectrometer  has  aeen  assembled  for  tne 
extended  free  spectral  range  up  to  100  cm  L  to  study  low  frequency  modes  in 
highly  transparent  solids.  The  systems  are  under  coraputor  control  using  a  CAMAC 
interface.  Aside  from  studying  inelastic  Brillouin  and  Raman  regimes  using  the 
interferometer,  the  electron!!,  system  enables  us  to  also  studv  quasi-elastic 
(Rayleigh)  scattering  by  digital-photon  counting  and  self-beating  techniques. 

Current  technological  interests  continue  to  motivate  a  search  for  candidate 
materials  which  exhibit  low  optical  losses  at  specific  wavelengths  for  utili¬ 
zation  in  fiber  optic  and  other  light  guiding  applications.  Considerations  have 
been  given  to  the  ultimate  theoretical  intrinsic-loss  mechanisms  and  the  prac¬ 
tical  extrinsic  losses  due  to  impurities  and  imperfections  in  crvstals  and 
glassy  materials.  The  ultimate  transparence  of  a  material  should  lie  in  a  valley 
between  the  intrinsic  electronic  absorption  edge  and  the  tail  of  the  muitiphonon 
transitions.  In  this  spectral  region,  the  ultimate  losses  should  be  due  to 
Rayleigh  and  Brillouin  scattering.  The  Brillouin  scattering  results  from  acousti 
cal  phonons  and  is  an  intrinsic  property  of  a  given  material.  The  Rayleich 
scattering  cf  light  results  from  nonpropagating  fluctuations  in  the  dielectric 
constant.  These  fluctuations  can  result  from  frozen-in  random  variations  in  the 
dielectric  constant  of  a  disordered  solid  as  well  as  fluctuations  resulting  from 
local  variations  in  strain  in  crystalline  solids;  such  extrinsic  scattering  will 
vary  from  sample  to  sample;  however  entropv  fluctuations  which  are  intrinsic  pro¬ 
perties  can  still  contribute  to  the  Rayleigh  scattering. 

Using  our  Rayleigh-Briilouin  scattering  system,  we  have  completed  an  ini¬ 
tial  study  of  scattering  from  single  and  polycrystalline  materials  in  coopera¬ 
tion  with  a  group  at  the  Hughes  Research  Laboratory,  to  probe  the  scattering 
mechanisms  in  highly  transparent  solids.  We  have  found  that  polycrystalline 
materials  scatter  more  strongly  than  single  crystal  materials.  However,  we  have 
not  as  yet  explicitly  associated  a  particular  elastic  scattering  process  which 
contributes  to  the  total  scattering.  Grain  boundaries,  dislocations  and  the 
decoration  of  these  structures  by  impurity-pinning  are  possible  candidates.  In 
addition,  our  results  indicate  that  even  in  pure  single  crystal  'KC-  samples,  ther 
is  more  scattering  than  predicted  for  an  ideal  crystal. 
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A  paper  entitled:  "Scattering  Losses  in  Single  and  Polycrystalline  Materials 
for  Infrared  Fiber  Applications"  is  to  be  published  in  the  volume  entitled  "The 
Physics  of  Optical  Fibers,  1980"  as  the  Proceedings  of  the  Am.  Ceramic  Societv 
Annual  Meeting,  Chicago,  April  1980.  A  copy  of  this  manuscript  is  enclosed  in 
the  Appendix. 


Fluoride  Glasses 

The  search  for  materials  which  exhibit  low  optical  losses  has  been  directed 
progressively  to  larger  wavelength  for  the  following  criterion.  Consideration  of 
the  various  contributions  to  intrinsic  scattering  indicate  that  the  minimum  opti¬ 
cal  loss  should  occur  at  the  intersection  of  the  Rayleigh  scattering  curve  on  the 
high  frequency  end  and  at  the  multiphonon  edge  at  low  frequency.  Since  the  Rayleigh 
scattering  losses  exhibit  a  A  wavelength  dependence,  it  is  clear  that  it  is 
desirable  to  operate  at  longer  wavelengths. 

It  has  recently  been  discovered  that  it  is  possible  to  prepare  a  series 
of  non-oxide  glasses  based  on  a  mixture  of  ZF,  and  H£F,  with  other  metallic  fluo- 
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rides.  In  the  conventional  preparation  of  these  glasses  residual  impurities  such 
as  OH  and  0  can  contribute  extrinsic  absorption  and  scattering  centers.  Because 
of  our  programmatic  interest  in  highly  transparent  materials,  we  have  turned  our 
attention  to  the  fabrication  and  optical  measurements  on  these  glasses.  In  collab¬ 
oration  with  a  group  at  the  Hughes  Research  Laboratory,  glasses  consisting  solely 
of  high  purity  ZnF, ,  ThF. ,  and  BaF„  have  been  synthesized  using  reactive  atmos- 
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phere  processing  (RAP)  techniques.  The  utilization  of  RAP  in  the  synthesis  cf 
these  glasses  has  shown  that  infrared  transparency  and  the  mechanical  strength  of 
these  halide  glasses  are  greatly  improved  by  eliminating  anionic  impurities  such 
as  OH  and  0  resulting  in  glasses  which  are  continuously  transparent  from  0.3  to 
7  Jjm  and  are  water  insoluable  and  unusually  strong  and  hard. 

This  work  has  been  published  in  "Infrared-Transparent  Glasses  Derived  from 
the  Fluorides,  Thorium  and  Barium,"  M.  Robinson,  R.  C.  Pastor,  R.  R.  Turk,  D.  P. 
Devor,  M.  Braunstein  and  R.  Braunstein,  Mat.  Res.  Bull.  L5,  735  (1980). 


Rayleigh-Brillouin  and  Raman  Scattering  in  Fluoride  Glasses 


In  order  to  evaluate  the  scattering  losses  in  these  glasses,  Ravleigh- 
Brillouin  and  Raman  measurements  were  performed.  By  measuring  the  ratio  of 
the  intensity  of  the  Rayleigh  to  Brillouin  scattered  light  -  the  so-called 
Landau  Placzek  ratio,  one  can  obtain  the  attenuation  due  to  scattering  in 
terms  of  calculable  scattering  losses.  Such  measurements  were  performed  on 
two  glass  systems:  HfF^-BaF?  -  ThF^  and  ZrF^-ThF^-BaF^ .  Preliminary  measure¬ 
ments  indicate  that  the  scattering  is  larger  for  the  HfF.  glass  than  the  ThF. 
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glass.  Whether  this  is  fundamentally  due  to  the  "fictive  temperature"  of 
these  glasses  or  is  extrinsic,  remains  to  be  seen  after  further  measurements 
on  a  number  of  samples  are  performed. 

In  addition,  Raman  scattering  measurements  were  made  on  these  systems  as  a 
means  of  identifying  the  optical  phonons  in  the  glasses  compared  to  the  phonons 
of  the  pure  constituents  of  the  glasses. 


Surface  Polaritons 

Optical  attenuated  total  internal  reflection  measurements  were  performed  on 
the  surface  polariton  of  n+  Silicon  and  Copper.  These  investigations  enable 
the  optical  contants  of  the  materials  to  be  determined  in  a  spectral  regime 
where  they  are  difficult  to  measure  by  conventional  means  and  in  addition,  the 
dispersion  relationship  of  the  surface  excitations  were  obtained.  The  general 

dispersion  relationship  for  surface  plasmon  at  an  air-solid  interface  is 

2  ,  , 

k-^  =  (oj'/c“)  (e(6j)/(l+e(w))  where  e(w)  is  the  complex  bulk  dielectric  func¬ 

tion  of  the  solid.  In  general,  there  will  be  two  physical  solutions  to  the  dis¬ 
persion  relationship  that  is,  one  where  the  frequency  is  real  and  k^  complex; 
and  the  other  where  (o  is  complex  and  k-Q  real.  Up  to  the  present  work,  there 
were  wide  discrepancies  between  theoretical  and  observed  dispersion  curves  ob¬ 
tained  by  attenuated  total  internal  reflection  techniques.  Previously,  investi¬ 
gators  had  attempted  to  interpret  their  data  using  real  values  of  U).  Hue  to  the 
apparent  discrepancies  between  theoretical  predictions  and  experimental  evalua¬ 
tions  in  the  real  frequency  phase,  we  examined  the  surface  polariton  dispersion 
relationship  for  various  models  in  the  complex  frequency  plane.  An  original 
method  was  developed  to  evaluate  the  optical  constants,  dielectric  functions 
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and  the  Fresnel  reflection  coefficients  by  redefining  them  as  functions  of 
complex  frequency.  By  evaluating  the  attenuated  total  integral  reflection 
dispersion  curves  in  the  complex  plane,  we  were  able  to  obtain  agreement  be¬ 
tween  theory  and  experiment  for  the  materials  measured  in  this  work.  Before 
this  analyses,  it  was  believed  that  the  backbending  observed  in  some  disper¬ 
sion  curves  would  occur  only  for  complex  wavevector  solutions  of  the  dispersion 
relationship.  Our  analysis  comparing  the  complex-4)  and  complex-k  solutions  of 
the  dispersion  relationship  showed  that  the  difference  between  the  two  calcu¬ 
lations  is  small. 

As  a  consequence  of  this  work,  we  believe  that  all  studies  of  surface  polar- 
itons  should  be  conducted  in  the  complex  frequency  plane  when  either  the  damp¬ 
ing  is  large  or  there  exists  interband  transitions  in  the  spectral  legion 
studied. 

The  above  work  has  been  completed  and  the  following  manuscripts  have  been 
submitted  for  publication: 

"Optical  Properties  of  the  Surface  Plasmon  of  n  Silicon"  -  P.  F.  Robusto 
and  R.  Braunstein  to  appear  in  Phys.  Stat.  Sol.  (b) ,  "Optical  Measurements 
of  Surface  Plasmons  in  Copper"  -  P.  F.  Robusto  and  R.  Braunstein,  to  appear 
in  Phys.  Stat.  Sol.  (b) . 


Publications  For  This  Period 


"Infrared-Transparent  Glasses  Derived  From  the  Fluorides  of  Zirconium, 

Thorium  and  Barium,"  with  M.  Robinson,  R.  C.  Pastor,  R.  R.  Turk,  D.  P.  Devor 
and  M.  Braunstein  (Mat.  Res.  Bull.  L5,  735  (1980). 

"Scattering  Losses  in  Single  and  Polycrystalline  Infrared  Materials  for 
Infrared  Fiber  Applications,"  with  J.  A.  Harrington,  M.  Braunstein  and  B.Bobbs 
(to  appear  in  Proceedings  of  Am.  Ceramic  Society  Annual  Meeting,  April  1980, 
to  be  published  as  Volume  on  Physics  of  Optical  Fibers,  1980). 


"Infrared  Wavelength  Modulation  Spectroscopy  of  Laser  Window  Materials,"  with 
R.K.  Kim  and  M.  Braunstein;  Laser  Damage  in  Optical  Materials;  1980,  Oct.  1980 
MBS  Boulder  Co.  ed  by  H.  E.  Bennett,  A.  J.  Glass  and  A.  H.  Guenther  to  appear 
in  MBS  Special  Publication. 

"Infrared  Transparent  Glasses  Derived  from  Hafnium  Fluoride,"  with  M.  Robinson 
R.  R.  Turk,  D.  P.  Devor  and  M.  Braunstein;  To  appear  in  Proceedings  of  Los 
Angeles  Technical  Symposium  Feb.  9-13,  1981,  North  Hollywood,  California,  the 
Society  of  Photo-Optical  Instrumentation  Engineers. 
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"Optical  Measurements  of  Surface  Plasmons  in  Copper,"  with  P.  Robusto  to  appear 
in  Phvs.  Stat.  Sol.  (b) 

"Optical  Properties  of  the  Surface  Plasroon  of  n+  Silicon,"  with  P.  Robusto  to 
appear  in  Phys.  Stat.  Sol.  (b) . 
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Scientific  Interactions 

We  have  continued  active  relationships  with  other  groups  in  the  field. 

Liason  was  maintained  with  the  synthesis  groups  at  the  Hughes  Research  Labora¬ 
tory  regarding  the  evaluation  of  alkali  halide,  alkali  earths  and  halide 
glasses  prepared  by  reactive  processing.  The  work  on  Hg^^Cd^Te  oxides  has 
benefitted  from  samples  and  technological  exchanges  with  the  Santa  Barbara 
Research  Center.  The  Si  on  sapphire  and  ion  implanted  samples  were  obtained 
from  the  Industrial  Products  Division  of  Hughes  at  Carlsbad.  The  work  on  the 
fluoride  glasses  is  a  cooperative  endeavour  with  the  Hughes  Research  Laboratory. 
In  addition,  discussions  on  this  project  have  been  held  with  3.  Bendow  of  RADC. 
We  have  had  extensive  technical  exchanges  regarding  shallow  doner  and  acceptor 
levels  in  molecular  beam  epitaxial  material  with  C.  W.  Litton  of  the  Air  Force 
Avionics  Laboratory  and  through  an  extensive  visit  by  Dr.  K.  Baja  to  UCLA.  In 
addition,  we  have  participated  in  the  ARPA  review  on  nucleation  of  crystal 
growth  and  the  properties  of  HgCdTe.  We  have  evolved  a  detailed  program  with 
the  Morris  Braunstein  Molecular  Beam  Epitaxial  group  at  the  Hughes  Research 
Laboratory.  This  program  entails  a  unified  correlation  with  the  Hughes  syn¬ 
thesis  activity  and  the  UCLA  study  of  shallow  and  deep  levels  in  MBE  lavers  of 
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INFRARED  WAVELENGTH  MODULATION  SPECTROSCOPY  OF  LASER  WINDOW  MATERIALS* 


RuDin  Sraunstein,  Ryu  <.  Kim 
Department  of  Physics 
University  of  California 
Los  Angeles.  California  50024 
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Morris  Braunstein 
Muqnes  Researcn  Laboratories 
Malibu.  California  90265 


We  have  previously  shown  that  the  infrared  wavelength  modulation  system  that  we  have 
developed  is  capable  of  measuring  the  complete  spectral  distribution  of  the  extrinsic 
absorption  in  nignly  transparent  solids  at  levels  of  10"*  cm*'  in  the  spectral  range  from 
2.5  to  12  microns  [)].'  The  previous  measurement  techniques  that  were  employed  for  KBr  and 
KCi  have  been  extended  to  other  crystalline  systems  wnicn  are  of  interest  in  light  guiding 
applications,  namely:  CaF2«  LiF,  NaCi,  NaF,  LaFg,  BaF2,  MgFj,  SrF2,  and  MgQ.  Rich  anq 
■  varied  absorption  structures  were  observed  in  all  of  these  crystals  enabling  an  identifica¬ 
tion  of  volume  and  surface  absorption.  Similar  dominant  bancs  are  observed  in  many  of  these 
substances  Indicating  the  presence  of  common  impurities  regardless  of  the  crystal  and  the 
origin  of  its  growth.  However,  varied  fine  structures  are  observed  in  different  crystals 
wnich  are  indicative  of  the  individual  characteristics  Of  the  chemistry  of  the  crystal 
preparation.  These  measurements  were  performed  in  laboratory  and  dry  N2  ambients  and 
readily  show  the  physisorption  and  desorption  of  surface  contaminants.  In  appropriate 
materials,  measurements  were  extended  into  the  intrinsic  multiphonon  region. 

Key  words:  Extrinsic  and  intrinsic  infrared  absorption;  volume  and  surface  impurities;  infrared  wave¬ 
length  modulation;  CaF^;  L* F ;  NaCi;  NaF;  LaF^;  BaF,;  MgF.,,  SrFj,  MgO;  laser  windows. 

1 .  Introduction 

The  requirements  for  high  power  laser  and  light  guiding  systems  continue  to  generate  a  need  for 
extremely  low  absorption  levels  in  infrared  transmitting  materials.  These  interests  continue  to 
stimulate  efforts  at  material  growth  and  surface  preparation  techniques  to  produce  materials  that 
approach  intrinsic  absorption  levels  in  appropriate  spectra)  regions.  In  practice,  in  the  highly  trans¬ 
parent  regions,  absorption  coefficients  are  observed  to  vary  greatly  from  sample  to  sample  indicating 
the  extrinsic  nature  of  the  absorption  processes.  To  assess  progress  in  improvements  in  material 
preparation,  sensitive  techniques  are  required  to  measure  low  levels  of  adsorption. 

To  determine  the  sources  and  nature  of  the  extrinsic  light  absorptions,  we  have  emoloyed  the 
technique  of  infrared  wavelength  modulation  in  the  spectral  region  from  2.5  to  12.0  urn.  we  nave  ore- 
viously  employed  this  technique  [1]  in  a  study  of  KBr  and  <C l  orepa-eo  by  different  crystal  growth 
techniques  and  were  able  to  identify  possible  surface  and  volume  impurities.  In  the  present  wort,  we 
have  extended  these  studies  to  survey  a  number  of  ootical  materials  commonly  used  in  laser  window 
applications,  namely:  CaF 2,  LiF,  NaCi,  NaF.  LaF3,  BaFj,  M9F2,  SrF2,  and  MgO.  The  continuous  spectral 
distribution  of  the  absorption  obtainable  by  wavelength  moduUiion  tecnniaues  at  'evels  of  10*5  cm-' 
facilitates  identification  of  volume  and  surface  absorption  tenters.  Althouqn  rich  and  varied  absorp¬ 
tion  structures  are  observed  in  all  of  these  crystals,  similar  dominant  bands  are  observed  in  many  of 
these  substances  indicating  the  presence  of  common  impurities  regardless  of  the  crystal  origin  ana  the 
concomitant  growth  techniques  used  to  prepare  the  materials. 

2.  Experimental  Techniques 

The  Infrared  wavelength  modulation  spectrometer  was  previously  described  J'].  It  consists  of  1 
system  for  sweeping  the  output  wavelength  of  a  monochromator  across  the  exit  slit  by  a  vibrating  out¬ 
put  diagonal  mirror  whose  depth  of  modulation  can  Be  continuously  controlled.  Careful  consideration 
is  given  to  the  subtraction  of  the  background  of  the  transfer  function  of  the  monochromator  wnicn  can 
include  atmosoheric  absorption,  the  soectral  response  of  the  optics  as  well  as  the  detector.  *0 
accomplish  this,  a  sample-tn-sample-out  procedure  is  used  to  obtain  the  background  intensity  without 
the  sample,  the  derivative  of  this  background,  the  transmitted  intensity  with  the  sample  present  mo 
its  derivative.  The  above  operations  are  performed  under  microprocessor  control,  'he  computation  of 
the  energy  derivative  of  the  absorption  as  well  as  Its  integration  to  display  absorption  peaks  are 
performed  on-line  using  a  P DP  11/20  computer.  Although  the  primary  data  in  the  fonn  of  the  energy 
derivative  of  the  sample  absorption  is  adequate  to  reveal  absorption  structures,  as  an  aid  for  com¬ 
parison  with  the  conventional  literature,  the  Integrated  results  are  reported  in  this  paper. 

3.  Experimental  Results  and  Analysis 

The  crystals  used  in  this  study  are  indicated  in  table  1,  tooether  with  suppliers  of  the  materials. 
Standard  procedures  were  used  to  prepare  the  surfaces  of  the  samples,  these  consisted  of  surface 
grinding,  polishing  with  a  slurry  of  Linde  A  in  ethanol  or  propanol,  ano  then  drying  under  a  neat 
iamo.  Prior  to  the  ootlca'  measurements  the  samples  were  stored  m  a  vacuum  desiccator.  The  optical 
measurements  were  performed  with  the  sawoles  in  the  ’aboratory  ambient  as  well  as  in  a  iry  N;  atmos- 
onere.  'he  data  for  the  various  materials  studied  are  shown  in  'igures  1.  3,  4,  S,  6,  S,  3,  "3 


'Suooorted  By  the  Air  ~orce  Office  of  Scientific  Researcn. 

'Figures  : n  oracKets  indicate  i’terature  references  at  the  end  of  this  paper. 
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The  laboratory  ambient  and  dry  No  atmosphere  runs  for  a  given  samcle  are  displayed  on  the  same  page  to 
facilitate  comoarison.  The  absorption  coefficients  at  the  upoer  'eft  hand  of  the  "'cures  are  the 
beginning  of  the  intrinsic  multipnonon  tails  and  are  shown  wnen  they  fall  within  the  spectral  range  of 
the  spectrometer. 

When  the  absorption  coefficient  is  sufficiently  hign,  as  in  the  multipnonon  region  it  is  possible 
to  measure  the  transmission  directly  as  well  as  its  derivative.  In  this  case,  it  is  possible  to 
determine  the  constant  of  integration  wnicn  is  necessary  to  obtain  the  absolute  value  of  the  absorp¬ 
tion  wnen  integrating  derivative  data;  tms  constant  can  be  used  to  calculate  the  absorotion  coeffi¬ 
cient  in  the  extrinsic  region  well  below  tne  multipnonon  tail.  The  constant  of  integration  can  also 
be  detained  from  a  laser  calorimetric  measurement  at  a  given  freouency.  The  integration  is  performed 
by  a  standard  trapezoidal  integration;  in  this  computation  the  depth  of  modulation,  wnich  in  tne  pres¬ 
ent  work  is  10  cm*',  is  inserted  as  a  parameter.  The  integrated  derivative  data  without  the  constant 
of  integration  is  shown  to  the  right  of  the  figures.  The  zero  of  tic  is  obtained  by  a  computer  cor¬ 
relation  procedure  to  deconvolute  the  background  derivative  [2],  The  zero  of  the  IK  wavelength  modu¬ 
lation  result  is  subsequently  normalized  to  an  absorption  coefficient  determined  by  laser  calorimetry 
shown  in  eacn  figure,  wnicn  effectively  supplies  the  constant  of  integration.  To  read  tne  actual 
absorption  coefficient  at  a  given  wavelength,  one  merely  adds  or  subtracts  the  value  of  AK  at  that 
wavelength  to  the  calorimetlve  value. 

The  absolute  value  of  the  absorotion  coefficient  can  be  slightly  ambiguous  in  providing  the  base¬ 
line  data  using  multi-mode  laser  calorimetry  since  the  mode  structure  of  tne  laser  can  straddle  some 
line  features  in  tne  absorption  spectrum  of  the  sample.  Obtaining  the  constant  of  integration,  i.e., 
the  baseline  correction  by  a  direct  absorption  measurement  in  the  same  wavelength  modulation  apparatus 
would  be  a  more  direct  procedure  as  when  tne  derivative  measurements  span  the  spectral  range  into  the 
multipnonon  region.  This  was  not  done  in  this  work  since  this  procedure  would  be  more  precise  if  the 
direct  aosorption  measurement  were  performed  above  and  below  the  region  where  the  baseline  is  to  be 
obtained.  However,  since  the  main  thrust  of  the  present  work  is  to  reveal  the  spectral  features  of 
the  absorption,  the  laser  calorimetric  normal ization  was  used  and  the  values  of  AK  from  the  integrated 
wavelength  modulation  data  displays  tne  fine  structure  excursions  above  and  below  the  calorimetric 
point. 

Figure  la  and  figure  lb  show  the  absorption  spectra  of  BaFj  taken  in  the  laboratory  ambient  and 
in  a  dry  ho  atmosphere,  resoectively.  The  richness  of  the  spectra  should  be  noted;  previous  laser 
calorimetric  measurements  using  a  few  discrete  lines  snow  only  a  broad  peak  between  1-5  jn  peaking  at 
3.4  jin.  Comparing  the  spectra  in  figures  la  and  lb;  it  should  be  noted  that  the  absorotion  within  the 
spectral  region  from  1600  to  4000  cm"'  decreases  after  the  sample  has  been  initially  in  the  laboratory 
atmosphere  and  is  subsequently  placed  in  a  dry  ho  atmosphere  for  two  hours.  Hhen  the  samole  is 
returned  to  the  laboratory  ambient  for  several  nours  the  original  structure  In  figure  la  is  recovered. 
The  fact  that  the  magnitude  and  the  linewidths  of  some  of  the  bands  cnange  in  a  dry  No  atmosphere 
indicates  that  a  portion  of  the  absorption  is  due  to  surface  onysisorotion.  The  structures  that  per¬ 
sist  even  wnen  tne  samoles  are  flushed  in  dry  N?  may  be  due  to  surface  chemisorption  species  or  volume 
'Impurities.  The  absorption  edge  starting  at  1050  car'1  and  Increasing  to  lower  energy  is  the  multi¬ 
phonon  absorption  tail. 

Crystals  containing  substitutional  molecular  anion  impurities  typically  exhibit  a  number  of  sharp 
infrared  aosorption  bands  corresponding  to  various  internal  modes  of  polyatomic  ions.  There  may  be 
small  shifts  in  frequency  and  increase  in  linewidths  depending  an  the  particular  host;  however, 
because  of  the  relatively  small  interaction  with  the  crystal  host,  the  vibrational  frequencies  can  be 
used  to  identify  a  particular  chemical  impurity  and  the  nature  of  its  bonding.  In  addition  to  bulk 
absorption  there  Is  the  possibility  of  surface  absorption  wmch  can  be  comparable  to  bulk  absorotion 
which  can  be  due  to  crystal  growth  techniques,  mechanical  or  chemical  polishes  and  chemical  cleaners. 

In  addition,  exposure  to  the  atmosphere  can  lead  to  a  deposition  of  a  conglomerat'on  of  impurities. 

Correlation-type  charts  have  been  published  wmch  tabulate  molecular-ion  vibrational  frequencies 
as  a  means  of  identification  of  possible  surface  and  volume  ‘mpurules  m  laser  windows  [3].  If  one 
considers  a  complete  frequency  overlay  of  all  possible  impurities  that  can  be  present  in  concentra¬ 
tions  of  0.1  opm,  one  would  expect  a  duasi-continuian  aosorption  throughout  the  2.5  to  '2.3  :.m  region 
due  to  the  overlap  df  the  lorentzian  tails  of  the  various  aosorption  bands.  However,  the  distinct 
absorption  bands  observed  in  the  present  work  indicates  the  presence  of  dominant  impurities.  Conse¬ 
quently,  rather  than  dwelling  an  the  multiplicity  of  possible  impurities,  we  shall  consider  tne  common 
molecular  anion  impurities  suen  as  the  metaborates ,  nitrides,  mydrocarbons ,  carbonates,  and  OH"  ‘ons 
that  are  revealed  by  conventional  absorption  spectroscopy  of  Pure  commercial  crystals  since  't  >s 
reasonaole  to  expect  that  they  may  still  oe  present  at  lower  levels  in  purified  materials. 

The  vibrational  frequencies  of  common  bonds  of  some  polyatomic  anions  are  shown  in  fiqures  2a  and 
2b  which  display  the  OH",  C-h  and  carbonate  groups  and  tne  metaoorates  and  nitrides,  respectively; 
figure  2a  will  be  used  primarily  In  the  following  analysis  since  these  impurities  seem  to  be  dominant; 
however,  bands  in  figure  2b  may  also  be  present. 

If  we  examine  figure  la  for  the  laboratory  atmosphere  data  In  conjunction  with  figures  2a  and  2b, 
it  is  possible  to  Identify  some  possible  impurities.  4  prominent  band  at  4.2  jii  can  be  associated 
with  the  ohysisorption  of  atmospheric  COj.  The  bands  in  the  4  to  6  jm  and  2.3  -m  -egions  can  be  lue 
to  atmospheric  water,  while  the  bands  in  the  3  to  4  jm  region  can  be  assnned  to  an  overlap  of  hydro¬ 
carbon  bands.  The  bands  ’n  the  6  to  3  jn  region  seem  to  be  associated  with  carbonates.  rhe  ’ act  *hat 
the  oands  in  the  4  to  6,  2.3  jm  and  3  to  4  jn  regions  are  'argely  due  to  pnysisorbed  species  's  con¬ 
firmed  by  the  decrease  in  absorption  in  the  dry  Nn  atmosonere  data  shown  in  fioure  'b. 

The  region  between  6  and  3  jii  in  BaF?  snows  considerable  structure  wnose  magnitude  varies  only 
slightly  with  dry  N->  'lusning.  4s  we  have  previously  indicated  by  the  examination  pf  f'gures  la  and 
lb  in  conjunction  with  ‘igures  2a  and  2b  these  structures  can  be  assoc'ated  w'th  tne  rhemisorptlon  of 
carbonates.  5hotoacoust“c  measuremerrs  on  3aFi  and  Srf?  using  a  CO  laser  that  was  tunable  to  discrete 
1‘nes  •»  the  6  to  3  urn  region  "4]  nave  revealed  step-1 'ke  structures  (n  the  surface  adsorption,  wirle 
our  measurwwnts  re*ej!  distinct  bands. 
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The  MgF2  laboratory  atmosphere  data  shown  in  figure  3a  reveal  a  cluster  of  bands  similar  to  those 
ooserved  in  8aF?.  Again  the  ubiquitous  CO?  at  A. 2  urn  is  evident.  The  2.3  and  4  to  6  um  regions 
reveal  the  'liquid"  water  and  OH-  bands,  while  the  3  to  1  region  reveals  the  possible  overlap  of 
C-H  bands.  The  6  to  3  bands  seen  in  3aF?  seem  to  be  absent  or  are  obscured  by  the  strong  water 

band  in  the  4  to  6  um  region.  The  dry  M?  data  in  figure  3b  snows  a  dramatic  suppression  of  ail  the 
above  bands.  The  remaining  structure  can  be  due  to  volume  or  chemisorbed  species;  the  rise  in  the 
neighborhood  of  2.3  and  4.3  can  ce  due  to  OH'  and  "liquid"  water,  respectively. 

The  SrFj  sample  in  laboratory  atmospnere  shown  in  figure  4a  is  very  similar  to  the  structure  seen 
in  SaFj  and  the  identifications  of  tne  bands  are  similar  to  the  3aF?  discussion.  However,  it  should 
be  noted  that  Che  dry  M2  atmosphere  run  in  figure  4b  snows  a  mucn  greater  decrease  in  the  structure 
compared  to  the  SaFj  dry  M2  run.  The  remaining  structures  in  figure  4b  can  oe  due  to  hulk  or  chemi¬ 
sorption  with  the  possible  identifications:  6  to  8  um  (carbonates),  3  to  4  _m  (C-H  bands),  2.3  „m 
(OH"). 

The  CaF2  laboratory  data  in  figure  5a  again  shows  a  similar  distribution  of  bands  as  exhibited  by 

the  3aF?  and  SaF2  sample  with  a  possible  similar  identifications  of  the  bands.  The  marked  decrease  in 

height  of  these  bands  in  the  dry  M2  is  again  evident  as  shown  in  figure  5b. 

The  distributien  of  the  bands  in  MaF  in  figure  6a  exhibits  the  ‘amillar  pattern  previously  dis¬ 
cussed  for  BaF2,  SrF2,  and  CaF2-  However,  figure  6b  snows  a  large  decrease  in  the  band  heights  in  tne 
2.5  to  4.0  um  region. 

The  LaFj  sample  in  the  laboratory  atmosphere  shown  in  figure  7b  reveals  the  4.2  um  CO2  band,  tne 
5  to  8  um  possible  carbonates  and  the  OH"  band  near  2.3  um.  However,  it  should  be  noted  that  the 
bands  in  the  3  to  4  um  region  due  to  C-H  vibrations  wnicn  have  been  prominent  in  3aF?,  MqFj,  SrF?, 

CaF2  and  MaF  seem  to  be  aosent  or  greatly  suppressed.  In  the  dry  M2  data  in  figure  7b  the  oeak  near 
2.8  urn  possibly  due  to  OH'  is  about  tne  same  height  as  in  figure  7a.  One  should  note  the  decrease  in 

the  overall  absorption  in  the  4  to  3  um  region  as  indicated  by  tne  large  negative  values  of  the  ZK. 

The  MaCZ  data  in  both  the  laboratory  and  dry  M?  atmosohere  data  in  figure  Sa  and  3b,  respectively, 
show  prominent  structures.  Although  tne  spectral  distribution  is  different  in  detail  in  both  of  these 
figures  one  can  still  identify  the:  6  to  12  um  (carbonates),  4.3  um  ("liquid"  water),  the  4.2  jn  ( CO2 > » 
3  to  4  um  (C-H)  and  the  2.3  um  (OH")  bands.  The  prominence  of  the  bands  in  both  ambients  is  consis¬ 
tent  with  the  greater  surface  activity  expected  for  MaCZ  as  compared  to  the  above  substances.  Although 
there  is  a  large  decrease  in  the  height  of  tne  bands  in  figure  3b.  as  shown  by  the  large  negative 
values  of  J.K,  most  of  tne  bands  are  still  evident  that  were  present  in  figure  3a. 

The  lIF  data  in  figures  9a  and  9b  snow  little  structure  indicating  that  LiF  is  the  least  surface 
active  of  all  the  substances  studied.  The  absence  of  tne  4.2  -m  232  band  wnicn  was  visible  ;n  all  the 
previous  suostarces  studied  including  <Br  and  <CZ  wnicn  was  previously  published  [1]  should  be  noted. 
The  only  band  wnich  seems  to  be  noticeaole  is  4.5  um  band  in  both  ambients  possibly  due  to  water. 

The  data  for  MgO  are  shown  in  figures  10a  and  10b.  The  most  prominent  features  seen  ;n  the  lab¬ 
oratory  ambient  are:  a  doublet  in  the  2.3  to  3.0  .»  region  (water  and  CM'S,  some  aspects  of  a  pos¬ 
sible  4.2  um  (CO2),  structure  in  the  4  to  6  urn  region  (liquid  H^O )  and  very  slight  structure  in  the 
3  to  4  um  region  (hydrocarbons).  It  is  interesting  to  note  that  of  the  doublet  in  the  laboratory 
ameient  it  is  only  tne  3.0  um  band  wnicn  survives  dry  M?  flushing  with  a  noticeaDle  suppression  of  the 
2.8  jn  band.  In  a  previous  study  of  MgO  [5]  impurity  bands  were  ODserved  between  tne  3.3  -m  and 
2.7  um.  Our  study  clearly  shows  that  the  3.0  -m  band  is  due  to  bulk  or  chemisorbed  OH"  wnile  the 
2.3  um  band  is  due  to  ohysi sorbed  OH". 

The  region  of  intrinsic  absorption  in  MgQ  shows  some  fluted  structures  in  the  multiphonon  aosoro- 
tion  tail.  The  rise  around  1600  cm*'  (6.3  um)  agrees  with  the  shoulder  previously  studied  the 
multipnonon  spectra  of  MgO  [6]  and  was  attributed  to  4T0  phonons  by  suitably  averaging  over  the  ois- 
oerson  curves. 


4.  Summary 

Figure  11  shows  a  "schematic"  representation  of  the  typical  spectra  seen  in  most  of  tne  samples 
in  the  laboratory  atmosphere;  this  "free-hand"  drawing  essentially  summarizes  the  dominant  bands  that 
are  observed!  An  examination  of  all  of  the  samples  studied  except  LiF  show  many  features  in  common 
despite  the  fact  that  a  cursory  glance  would  emphasize  the  differences  in  detail.  The  6  to  S  .« 
(carbonates),  the  4-6  -m  (liquid  water),  the  4.2  um  (CO2 ) ,  the  overlap  of  bands  in  the  3  to  4  ,»  region 
(hydrocaroons)  and  the  3  and  2.3  jn  bands  (OH"  bound  to  water  and  free  3H")  3re  seen  in  most  of  the 
samples.  In  general,  the  ohyslsorbed  surface  character  of  these  banos  Is  confirmed  by  their  marked 
suppression  in  a  dry  M2  atmosphere. 


5.  Conclusion 

He  have  measured  the  continuous  soectral  distribution  of  the  absorption  In  a  number  of  materials 
of  Interest  for  laser  window  applications  In  the  spectral  reqlon  from  2.5  to  12.9  _m.  This  *jrk  nas 
enabled  an  identification  of  volume  and  surface  contaminants  where  their  existence  could  only  be  con¬ 
jectured  from  previous  discrete  thermocouple  laser  calorimetry  and  onotoacoustic  calorimetry  measure¬ 
ments.  It  is  interesting  to  note  that  although  the  diverse  chemistry  and  crystal  jrowtn  techniques 
were  employed  in  the  preparation  of  these  materials,  similar  dominant  bands  are  observed  in  most 
samoles  which  are  due  to  ohysi sorption,  volume,  and  surface  chemisorption.  It  nas  been  possible  to 
identify  the  type  of  bonding  of  the  particular  constituents  present,  "he  surface  character  of  some  of 
the  apsorption  at  levels  of  10"1  to  1 0” 5  cm-i  indicates  the  ’mportance  of  aer«orming  measurements  m 
the  same  ambient  when  comparing  ’ow  level  absorption  measurements  by  various  sensitive  techniques. 

"he  sensitivity  bf  tne  infrared  wavelength  modulation  technique  ‘ndlcates  that  a  fraction  of  a  nono- 
'ayer  of  surface  adsorbed  species  can  be  detected.  Experiments  in  controlled  ambients  and  on  well 
characterized  surfaces  usinq  this  technique  can  be  rewarding  'or  the  study  of  the  onysics  and  chemistry 
of  surfaces,  "he  dramatic  changes  ‘n  the  physisorbed  spectra  in  tne  substances  studied  ’noicate  that 
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a  continuous  flow  of  dry  along  tne  outer  surfaces  of  optics  may  nave  some  rewards  in  high  power 
laser  systems.  ‘ 
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Tadie  1.  Origins  of  tne  samoles  used  in  this  work. 
3a f2  -  HARSHAW 
MgF,  -  OPTOVAC 

Srf,  -  Hughes  Research  Laboratories  -  Press-forged 
from  RAP  Boule 

Ca F2  -  HARSHAW 

MaF  -  QPTQVAC 

LaFj  -  Hughes  Research  Laboratories  -  rie*HF  cast 
NaCl  -  Hughes  Research  Laboratories 
l if  -MEU.ER 


MgO  -  OPTOVAC 


FIGURE  CAPTIONS 


Figure  la.  Wavelength  modulation  spectra  of  3aF2 
in  the  laboratory  atmospnere. 

Figure  lb.  Wavelength  modulation  spectra  of  3aF? 
in  a  dry  atmosonere. 

Figure  2a.  Correlation  chart  of  carbonates, 
C'H,  0-H  and  v<ater  frequencies. 

Figure  2b.  Correlation  chart  of  metaborate  and 
nitride  frequencies. 

Figure  3a.  Wavelength  modulation  spectra 
of  MgFj  in  the  laboratory  atmosphere. 

Figure  3b.  Wavelength  modulation  spectra 
of  MgF,  in  a  dry  Nj  atmosphere. 

Figure  da.  Wavelength  modulation  spectra 
of  SrFj  in  the  laboratory  atmosphere. 

Figure  4b.  Wavelength  modulation  spectra 
of  SrF^  in  a  dry  N2  atmosphere. 

Figure  5a.  Wavelength  .modulation  spectra 
of  CaF^  in  the  laboratory  atmosphere. 

Figure  5b.  Wavelength  modulation  spectra 
of  CaF^  in  a  dry  atmospnere. 

Figure  6a.  Wavelength  modulation  spectra 
of  MaF2  in  the  laboratory  atmosphere. 

Figure  6b.  Wavelength  modulation  spectra 
of  NaF2  in  a  dry  fi^  atmospnere. 

Figure  7a.  Wavelength  modulation  spectra 
of  LaF^  in  the  laboratory  atmosphere. 

Figure  7b.  Wavelength  modulation  spectra 
of  LaF^  in  a  dry  N,  atmosphere. 

Figure  3a.  Wavelength  modulation  spectra 
of  NaCZ  in  the  laboratory  atmosphere. 

Figure  8b.  Wavelength  modulation  soectra 
of  NaCZ  in  a  dry  it,  atmosphere. 

Figure  9a.  Wavelength  modulation  soectra 
of  LiF  in  the  laboratory  atmospnere. 

Figure  9b.  Wavelength  modulation  soectra 
of  LiF  in  a  dry  N,  atmosphere. 

Figure  10a.  Wavelength  modulation  soectra 
of  MgO  in  the  laboratory  atmosphere. 

Figure  10b.  Wavelength  modulation  spectra 
of  MgO  in  a  dry  h,  atmospnere. 

Figure  11.  "Schematic''  representation  of 
the  typical  spectra  seen  In  most  samples 
in  the  laboratory  atmospnere. 
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ABSTRACT 

Polycrystalline  fiber  waveguides,  fabricated  from  infrared 
transparent  solids  such  as  KRS-5  and  KC1,  have  measured  losses 
much  greater  chan  conventional  silica  fibers.  One  major  source  of 
these  losses  is  scattering  from  grain  boundaries  present  in  the 
polycrystalline  fibers.  To  improve  the  optical  transmission  of  our 
infrared  waveguides,  we  have  studied  the  losses  due  to  scattering 
in  single  and  polycr^stalline  materials  which  are  suitable  for 
fabrication  into  infrared  transmissive  waveguides. 
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INTRODUCTION 


Optical  fiber  waveguides  made  from  crystalline  materials  such 

as  KRS-5  (TIBrI),  TIBr,  AgCl,  and  KC1  have  been  used  for  a  variety 
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of  10.6-ym,  CC?  laser  applications.  The  losses  in  the  current 

IR  fibers,  however,  are  high  (lowest  loss  measured  is  300  dB/km  at 

10.6  uni  in  KRS-5)  and  applications  in  sensor  and  laser  power 

delivery  systems  have  been  limited  to  snort.  (1-  to  2-m)  lengths  of 
2-4 

riber.  This  measured  IR  fiber  attenuation  is  considerably 

higher  than  that  predicted  theoretically  for  these  and  related  IR 
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transparent  materials.  ’  In  tact,  Gentile  et  al.  and  van  Uitert 
and  Wemple^  have  shown  that  these  materials  have  projected  losses 
as  low  a 10  ^  dB/km  near  5  ym.  To  develop  this  ultra-low- loss 
potential,  for  such  applications  as  long-distance  communication 
links,  requires  a  careful  analysis  of  the  nature  of  the  attenuation 
mechanisms  present  in  IR  transparent  waveguides.  In  this  paper,  we 
address  the  contribution  of  scattering  to  the  total  attenuation  in 
bulk  materials  that  have  the  potential  for  being  fabricated  into 
highly  transparent  fibers. 

The  attenuation  mechanisms  present  in  low-Loss  solids  are 
illustrated  in  Fig.  1  for  fused  silica.  At  the  shortest  wavelengths, 
electronic  processes  (Urbach  tail)  contribute  heavily  to  the  total 
loss.  At  the  IR  wavelengths  of  interest,  however,  two  mechanisms  — 
scattering  and  multiphonon  absorption  —  have  been  identified  as  the 
ultimate,  limiting  loss  processes. “  In  Fig.  1,  the  curves  for 


scattering  (which  decreases  as  A  with  increasing  wavelength)  and 
lattice  (multiphonon)  absorption  (which  increases  exponentially 
with  increasing  wavelength)  cross  to  yield  a  minimum  in  the  total 
attenuation.  For  fused  silica,  this  minimum,  which  is  about 
0.25  dB/km  at  1.6  um,  has  been  achieved  in  kilometer-long  fibers. ^ 
For  certain  crystalline  as  well  as  special  glassy  solids,  minima 
occur  near  5  um  with  projected  losses  well  below  the  intrinsic 
losses  measured  in  silica  (projected  losses  are  given  in  the  next 
section) . 

The  total  attenuation  coefficient  may  be  written  as  the 

sum 


■a 


T 
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(1) 


where  and  are  the  contributions  due  to  scattering  and  absorp¬ 
tion,  respectively.  Each  term  in  Eq.  1  can  be  measured  independently, 
thus  allowing  the  individual  mechanisms  contributing  to  the  overall 
optical  loss  in  solids  to  be  studied.  For  example,  otj,  can  be 
obtained  from  standard  spectroscopic  and  fiber  insertion  loss 
measurements  wnile  laser  calorimetry  has  been  used  verv  successfully 
to  determine  residual  absorption  ot  in  weakly  absorbing  materials. 

A 

The  scattering  terms  a,.  has  not  been  as  well  studied.  Measurements 
using  Integrating  spheres  for  both  bulk  and  fiber  materials  are 
generally  used  to  obtain  a  total  integrated  scattering  (TIS)  loss. 
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These  methods,  however,  have  the  disadvantages:  being  unable  to 
distinguish  among  the  various  individual  scattering  mechanisms  con¬ 
tributing  to  TIS.  To  elucidate  the  various  scattering  mechanisms 
as  well  as  to  obtain  a  value  for  oi  ,  we  have  chosen  to  study  the 
light  scattering  spectra  of  solids.  These  spectra  are  composed  of 
elastically  (Rayleigh)  scattered  light  that  results  from  various 
nonpropagating  fluctuations  in  the  materials  index  of  refraction 
and  inelasticallv  (3rillouin)  scattered  light  that  results  from 
the  interaction  of  light  and  the  thermal  motion  of  ions  (sound 
waves).  Although  Rayleigh-Brillouin  (RB)  spectra  have  been  used  to 
measure  i ^  in  glasses,^  this  technique  has  not  been  expressly  used 
before  to  study  scattering  losses  in  single-  and  poly-crystalline 
materials.  In  this  work,  we  have  measured  RB  scattering  at  90°  in 
bulk  single-  and  poly-crystalline  KC1.  As  discussed  in  the  next 
section,  we  expect  very  little  Rayleigh  scattering  in  single- 
crystal  materials;  for  polycrystalline  samples,  however,  intuition 
suggests  that  the  residual  strain  and  grain  boundaries  associated 
with  the  hoc  press-forged,  polycrystalline  material  should  lead  to 
larger  amounts  of  scattering.  Our  preliminary  results  supoort  this 
presumption,  but  we  have  not  yet  been  able  to  account  for  the 
source  of  each  scattering  mechanism  contributing  to  the  R3  spectra. 
GENERAL  CONSIDERATIONS 

The  Limiting  attenuation  mechanisms  in  transparent  solids  are 
scattering  and  auitiphonon  absorption.  These  contributions  nave 


been  considered  by  several  investigators  in  the  context  of  projecting 
future  ultra- low- loss  materials  for  the  next  generation  of  fiber 
waveguides.  Van  Uitert  and  Wemple^  have  studied  the  potential  of 
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ZnCl^  (a  glass  former)  while  Gentile  et  al.“  and  Pinnow  et  al. 
have  concentrated  on  the  crystalline  materials.  In  Fig.  2.  we  con¬ 
sider  the  projected  transmission  for  KRS-5  and  KC1  and  compare  these 
predictions  to  fused  silica.  The  curves  in  Fig.  2  show  the  charac¬ 
teristic  V  shape  resulting  from  the  crossing  of  the  scattering 
(short  wavelength)  and  multiphonon  (long  wavelength)  attenuation 

mechanisms.  The  only  scattering  mechanism  assumed  in  calculating 
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the  X  -dependent  scattering  curve  for  KRS-5  and  KC1  was  Brillouin 
scattering  (ot_  -  see  below).  As  mentioned  above,  the  V-shaoed 
curve  for  silica  has  essentially  been  traced  out  experimentally 
and  thus  silica  fiber  losses  are  now  intrinsic.  For  KC1  and  KRS-5 
(as  well  as  for  many  other  non-oxide  ionic  solids),  however.  Fig.  2 
shows  the  extremely  low  loss  potential  for  these  materials  near 
5  um. 

To  determine  from  RB  light  scattering  experiments  requires 
a  careful  measurement  of  the  intensity  of  both  the  Brillouin-  and 
Rayleigh-scattered  light.  Brillouin  scattering  results  from  light 
that  has  been  inelastically  scattered  (Bragg  scattering)  from 
acoustical  phonons  (sound  waves).  The  frequency  of  the  scattered 
light  is  Doppler  shifted  from  the  frequency  of  the  laser  light 
u.  by  an  amount  if. : 
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U"nV  sinCv/d) 


where  n  is  the  index  of  refraction  of  the  medium,  v  is  the  velocity 

of  sound,  ■?  is  the  scattering  angle,  and  V  is  the  vacuum  wavelength 

o 

of  the  lighc.  This  frequency  shift  has  been  well  scudied  in  alkali 
g 

halides.  The  Rayleigh-scattered  light  (central  maximum  at  a,  )  is 
due  to  scattering  of  light  from  nonpropagating  fluctuations  in  the 
dielectric  constant.  In  glasses,  mechanisms  which  give  rise  to 
these  fluctuations  include:  density  variations  resulting  from  the 
frozen-in,  random  variations  in  dielectric  constant  inherent  in  a 
disordered  solid;  concentration  fluctuations  resulting  from  the 
local  compositional  variations  present  in  mixtures;  and  entropy 
fluctuations  resulting  from  temperature  variations.7  Of  the  three, 
only  entropy  fluctuations,  which  are  very  weak,  would  be  present 
in  an  ideal  single  crystal. 

To  obtain  ig,  we  first  evaluate  the  intensity  ratio  of  the 
Ravleigh  (I  )  to  the  total  Brillouin  (2I„)  scattered  light.  This 

A  D 

ratio  is  called,  based  on  its  use  in  light  scattering  in  liquids, 
the  Landau-Placzek  ratio  R^  and  is  defined'  as 
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Strictlv  sneaking,  since  '  and  thus  R. „  ire  related  to  specific 
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scattering  mechanisms  <  such  is  those  discussed  above  for  glasses). 


measured  values  of  R, p  are  generally  regarded  as  a  property  of  a 

given  material  (e.g.,  fused  SiO.,  has  an  R^p  13,  while  33K,,0-b7S  i0_, 
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has  an  jv  10).  In  our  experiments,  we  aeasure  I  without,  in 
general,  knowing  the  specific  mechanisms  contributing  to  the 
Rayleigh  component  of  scattered  light.  Therefore,  we  should  more 
appropriated  speak  of  an  effective  Landau- P lac zek  ratio  with  I 

£\ 

representing  the  Intensity  of  the  central  maximum. 

The  measured  R^p  is  then  used  to  calculate  a^ ,  as  described 
10  7,  h 

by  Pinnow  et  al.  and  others,’  from  the  relationship 


■*s  =*  ^(R^p  +  I)  .  (3) 

where  is  the  small  residual  attenuation  coefficient  due  to 
Brillouin  scattering  alone.  It  is  given  by 


a 
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(I) 


where  is  the  photoelastic  (Pockels)  coefficient,  T  is  the 

temperature,  k^  is  Boltzmann's  constant,  and  3„  is  the  isothermal 

compressibility.  For  ideal  single  crystals,  I  2i  Chus  R^p  j  0 

and  (from  Eq.  3)  a  a,  a_.  This  leads  to  the  scattering  curves  in 
b  a 

Fig.  1  for  KR.S-5  and  XCi,  which  were  calculated  from  Eq .  1  alone, 
while  the  scattering  for  fused  siiic3  was  calculated  from  Eq .  3 

using  R^p  = 


23. 
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An  interesting  feature  of  measured  pS  for  crystalline  solids 
is  their  dependence  on  polarization.  For  our  single-crystal 
measurements  in  KC1,  a  [100]  crystal  orientation  was  used  for  most 
samples  studies.  With  incoming  light  along  the  [100]  direction 
and  scattered  light  along  the  [010],  the  polarization  of  the  inci¬ 
dent  and  analyzed  scattered  light  was  either  vertical  (V)  or 
horizontal  (H)  with  respect  to  the  scattering  plane.  For  this 
geometry,  we  may  determine  the  intensities  of  the  Brillouin  compo¬ 
nents  from  the  selection  rules  for  the  rock,  salt  structure  (0,  ) 

h 

and  the  differential  cross  section. ^  The  intensities  I  and 

allowed  vibrational  modes  for  the  various  polarizations  (phonon 

momentum  q  along  [110]  direction)  are  summarized  in  the  matrix 

given  in  Table  la.  Using  the  known  photoelastic  (?^9  and  J  and 
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elastic  (C^,  C^,,,  and  C^)  constants  for  KC1,  we  calculate  the 

matrix  elements  in  Table  la  and  give  them,  normalized  to  I  (the 

an 

weakest  intensity),  in  Table  lb.  Table  lb  shows  that  the  intensi¬ 
ties  of  the  longitudinal  modes  (frequency  equal  to  15.3  GHz)  in  W 
polarization  are  much  stronger  than  those  in  the  other  polariza¬ 
tions  (transverse  mode  frequency  equal  to  7.2  GHz).  This  means 
that  measured  Landau-Placzek  ratios  mav  vary  greatlv  depending  on 
crystal  orientation  and  polarization,  with  the  smallest  R^p  occur¬ 
ring  for  the  W  polarization. 


LIGHT  SCATTERING  MEASUREMENTS 

The  RB  spectra  were  recorded  using  as  a  source  an  Ar-ion  laser 
(Spectra  Physics  Model  165)  delivering  10  to  300  mW  of  single-line 
power  and  a  PZT-scanned  Fabry-Perot  spectrometer.  The  experimental 
set-up  is  shown  in  the  block  diagram  in  Fig.  3.  The  spectrometer 
is  a  Burleigh  Instruments,  Inc.,  actively  stabilized  Fabry-Perot 
interferometer  with  its  associated  photon  counting  electronics. 

Data  acquisition  is  provided  by  a  512-channel  Tracor-Northern 
multichannel  analyzer  (MCA).  This  MCA  has  proven  essential  for 
obtaining  good  S/N  ratios  for  these  crystals  (between  2000  and 
20,000  scans  are  generally  accumulated  for  each  spectrum) . 

To  improve  the  contrast  of  the  spectrometer,  we  added  a  three- 
pass  attachment  to  the  interferometer.  This  allowed  us  to  readily 
detect  the  Brillouin  components  in  the  polycry staliine  materials 
where  the  Rayleigh  scattering  is  mors  intense.  Fig.  4  shows  the 
results  of  a  90°  scattering  measurement  on  polycrystalline  KCi 
doped  with  1.75*  RbCl  taken  with  the  light  passing  once  (one-pass) 
or  thrice  (three-pass)  through  the  interferometer.  The  improved 
contrast  in  the  three-pass  case  is  obvious.  Note  in  particular 
the  resolution  of  the  transverse  modes  in  the  three-pass  case  tnat 
has  been  lost  in  the  Rayleigh  wing  in  the  one-pass  case.  Clearly, 
the  three-pass  arrangement,  which  has  been  used  in  all  our  measured 
data,  is  necessary  to  obtain  a  reliable  R,p.  The  resolution  of 
the  triple-pass  Fabry-Perot  spectrometer  was  0.03  cm  ^  (finesse 
equal  to  30). 
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RESULTS  OF  LIGHT  SCATTERING  EXPERIMENTS 

Measurements  of  scattering  losses  were  made  at  138.0  nm  in 

bulk  single  and  polycrystalline  SCI.  The  XC1  was  reactive  atmos- 

X  3 

phere  process  (RAP)  single-crystal  material  which  was  either 
used  in  oriented  single-crystal  (SC)  or  in  hot-press  for  zee, 
polycrystalline  (?)  form.  The  RB  scattering  data  shown  in  Pig.  5 
are  for  pure,  single-crystal  KC1  oriented  as  shewn  in  the  insert. 
From  che  data,  we  can  see  the  intense  stokes  (S)  and  anti-scokes 
(AS)  longitudinal  Brillouin  components  in  the  TV  polar izarion.  The 
Brillouin  components  are  also  seen  to  become  weaker  in  7K  < or  HV' 
and  HH  polarization,  which  is  in  qualitative  agreement  with  the 
results  stated  in  Table  1  for  this  scattering  geometry.  Similar 
results  were  obtained  for  other  RAP-grown  XCi  (SC)  although  the 
R.  pS  were  found  to  vary  somewhat  from  sample  to  sample. 

The  data  for  polycrystalline  XC1  (average  grain  size,  10  cm)  ar 
shown  for  two  polarizations  in  Fig.  6.  We  can  see  from  these  data 
the  intense  Rayleigh  scattering  typical  of  our  polvcrystalline 
samples.  In  this  sample,  we  also  note  the  presence  of  only  longi¬ 
tudinal  modes  in  W  polarization  and  transverse  modes  in  7H 
polarization.  This  leads  us  to  conclude  that  we  have,  by  chance, 
illuminated  an  axis  of  high  symmetry  in  this  sample.  Specif icallv , 
it  would  appear  from  che  data  that  the  crystallites  are  oriented 
along  the  [100]  direction  for  this  particular  experiment.  In 

to  observe  an  admixture  of  L  and  T  modes 
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general,  we  would  expec t 


consistent  with  a  random  oriencation  of  polycrysballine 
samples. 

Effective  Landau-Hiuczek  ratios  have  been  calculated  for  these 
two  samples  from  the  data  in  Figs.  5  and  6.  These  data  are  sum¬ 
marized  in  Table  2.  For  the  single-crystal  :<C1,  R^p  is  lowest  tor 
Che  W  polarization,  as  discussed  above,  but  the  R^pS  for  other 
polarizations  do  not  scale  with  the  predicted  Brillouin  intensi¬ 
ties  (see  Table  1).  This  is  due  to  the  polarization  dependence  of 
the  Rayleigh-scattered  light.  Our  measurements  of  I  indicate  that 
the  Rayleigh  scattered  light  is,  as  expected,  most  intense  for  the 
W  polarization.  For  this  polarization,  1^  is  approximately 
3  to  10  times  stronger  than  when  measured  under  HV,  VH,  or  HH  con¬ 
ditions.  More  detailed  studies  of  the  depolarization  ratio  might 
give  insight  into  the  nature  of  the  static  defect  contribution  to 
elascic  scattering.  In  general,  the  values  for  R,  „  obtained  for 
polycrystalline  KC1(P)  are  higher  than  those  for  KCl(SC)  (see 
Table  2).  Again,  the  W  polarization  for  this  unoriented  poly 
sample  yields  the  lowest  R,  and  the  HH  case  yields  the  highest 
R,  P  (this  trend  was  also  seen  in  anot.ier  KC1(?)  sample). 

Table  2  also  gives  the  values  of  calculated  from  Eqs.  3 
and  4,  which,  for  KC1  at  483.0  nm,  reduce  to, 

a  =  1.5  x  10  ^  (R,  -  1)  cm 
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These  attenuation  coefficients  due  to  scattering  may  be  compared 

to  the  absorption  coefficients  obtained  for  KC1  by  laser  calorimet- 

14 

ric  measurements  at  438.0  nm  by  Harrington  et  al.  They  found 

2.  i  3  x  10-4  cm-1  for  KC1(P)  and  a  ,<  2  x  lCf°  cm" L  for  KCl(SC). 

•»e  see,  therefore,  that  the  total  attenuation  X_(=  i  )  is 

T  iy  S 

largely  due  to  scattering. 

The  nature  of  the  Rayleigh  scattering  in  our  crystalline 
samples  is  not  completely  understood.  The  results  indicate  sub¬ 
stantial  elastic  scattering  beyond  that  predicted  above  theoretical! 
from  entropy  fluctuations.  One  may  consider  this  excess  parasitic 
scattering  as  arising  from  mechanical  and  chemical  defects  in  the 
crystal.  Although  all  the  KC1  has  been  RAP  purified,  we  cannot 
rule  out  different  amounts  of  chemical  impurities  in  each  sample. 

It  is  also  evident  that  residual  strain  is  present  in  both  single 
and  polycrystalline  samples  (observed  as  birefringence  in  crossed 
polarizers).  This  strain,  which  one  would  intuitively  expect  to 
be  greater  in  the  not-forged  polycrystalline  materials,  can  lead 
to  substantial  elastic  scattering.  In  particular,  one  would  sus¬ 
pect  grain  boundaries  as  a  potentially  strong  source  of  Rayleigh 
scattering  because  impurities,  voids,  and  high  strain  would  be 
more  prevalent  in  these  areas. 

In  our  poiycrystalline  samples,  it  was  not  possible  to  examine 
the  scattering  from  a  single  grain  boundary  because  the  average 
grain  sice  (10  urn)  was  much  smaller  than  the  scattering  volume 


(cylindrical  volume,  350  um  long  by  30  in  diameter) .  Brody 
et  al. , ^  however,  were  able  to  study  elastic  scattering  from  a 
single  grain  boundary  in  polycrystalline  calcium  fluoride 
(Irtran-3).  In  their  measurements,  the  grain  size  (150  um)  was 
larger  chan  the  scattering  volume  and  they  found  essentially  no 
difference  in  the  intensity  of  light  scattered  from  within  a 
crystallite  to  that  scattered  from  a  volume  containing  a  grain 
boundary.  One  might  conclude  from  their  results  that  grain 
boundaries  do  not  contribute  to  elastic  scattering;  however,  to 
assess  the  effect  of  grain  boundaries  properly,  one  must  study 
high-purity  material  with  a  low  Rayleigh  background.  Our  future 
experiments  will  look  more  closely  at  the  effect  as  well  as  the 
importance  of  residual  strain  on  the  Rayleigh  scattered  light. 
SUMMARY 

Our  RB  scattering  studies  on  bulk  KC1  have  been  an  initial 
attempt  to  probe  the  mechanisms  responsible  for  scattering  losses 
in  highly  transparent  materials.  We  have  found  that  polycrystal¬ 
line  materials  scatter  more  strongly  than  do  single-crystal 
materials,  but  we  have  not  as  yet  been  able  to  explicitly  associ¬ 
ate  a  particular  elastic  scattering  process  with  a  finite  con¬ 
tribution  to  the  total  scattering.  Our  results,  however,  do  indi¬ 
cate  that,  even  in  these  very  pure  KC1  samples,  there  is  more 
scattering  than  predicted  for  the  ideal  KC1  crystals  and  thus 
future  fiber  waveguides  from  these  low-loss  materials  may  be  limited 
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more  by  scattering  chan  absorption  losses.  Consequently,  it  is 
important  to  study  scattering  in  poly-  and  single-crystal  materials 
to  understand  the  origin  of  the  Rayleigh  scattering  and  thereby  to 
determine  methods  of  minimising  this  contribution  to  the  attenuation. 
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TABLE  I 


(a) 


(b) 


Intensities 

of  3rillouin 

Components  for  a 

Rocksalt  Crystal 

Mode 

1W 

IVH  »  IHV 

IHH 

(110) 

L 

4?‘i2 

0 

4P44 

CU  +  C12  +  C44 

C11  +  C12  +  C44 

(lib) 

T 

0 

0 

0 

(001) 

T 

0 

P44 

0 

2C44 

Mode 

IW 

IVH 

IHH 

(110)  L 

133 

0 

1 

(lIO)  T 

0 

0 

0 

(001)  T 

0 

4.8 

0 

TABLE  II 


Effective  R, _  and  Calculated  Values  of  a. 


Structure 

Polarization 

*LP 

(cm*1)  1 

Single 

W 

26.2 

4.1  x  10~5 

Crystal 

VH 

33.5 

5.9  x  10  * 

(SC) 

HV 

44.9 

6.9  x  10'5 

HH 

276.4 

4.2  x  10~4 

Polycrystalline 

W 

496 

7.5  x  10“4 

(P) 

VH 

330 

-3  « 

1.2  x  10  ' 

HV 

1598 

-  i 

2.4  x  10  1 

HH 

>5000 

>7.5  x  10"3  1 
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-  Schematic  representation  of  optical  attenuation  mechanisms 
in  solids. 

-  Projected  losses  for  some  IR  transparent  materials. 

4 

Short  wavelengths  are  bounded  by  l/\  -scattering  and 
long  wavelengths  by  multiphonon  absorption. 

-  Block  diagram  of  scanning  Fabry-Perot  interferometer  for 
measuring  RB  spectra. 

-  Improved  contrast  in  RB  spectra  of  poly  KC1  +  1.752  RbCl 

through  the  use  of  3-pass  attachment  on  Fabry-Perot 
Interferometer.  The  spectra  were  all  recorded  using  the 
3-pass  configuration.  , 

-  Scattering  spectra  for  single-crystal  KC1  for  three  dif¬ 
ferent  polarizations.  Landau-Placzek  ratios  are  given 
in  Table  2  for  these  data. 

-  Scattering  spectra  for  polycrvstalline  K.C1  (average  grain 
size  is  10  urn).  Landau-Placzek  ratios  are  given  in 
Table  2  for  these  data. 
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ABSTRACT 

Glasses  consisting  solely  of  high-purity  ZrF^.ThF^,  and  3a?2  have  been 
synthesized  using  reactive  atmosphere  processing  (RA?)  techniques. 

RAF  of  the  individual  components  and  the  molten  material  with 
anhydrous  HF  and  CC14  is  described.  The  glass  molds  easily  at  312°C 
and  1920  psl  with  a  high-fidelity  replication  of  the  die  surface. 

The  glass  is  water  insoluble,  unusually  hard  ana  strong,  and  con¬ 
tinuously  transparent  from  0.3  to  7  „m. 


Introduction 

During  the  past  several  years,  considerable  effort  has  been  airected  to  the 
preparation  of  exceptionally  hign  purity  halide  materials  for  a  multitude  of 
applications  ranging  from  laser  windows  to  fiber  optics.  RAP  has  been  suc¬ 
cessfully  developed  for  crystalline  XC1  to  bring  about  the  reductions  the 
10.6-um  absorption  coefficient  from  cm-^-  (Ref.  I)  to  ilO-14  cm'^  (Ref.  2). 

RAP  methods  increased  the  rupture  strength  by  one  order  of  magnitude  and 
greatly  reduced  surface  attack  by  ambient  humidity. 

RAP,  when  applied  to  the  single-crystal  growth  of  rare  earth  fluorides  i3), 
totally  eliminated  many  solid-solid  phase  transitions  thought  to  be  intrinsic 
to  these  compounds.  The  transitions  are  triggered  by  low  concentrations  of 
0H~;  minimizing  the  OH-  level  by  atmosphere  control  permitted  large  single 
crystals  (>1  cm)  to  be  grown  from  melts  of  fluorides  previously  thought  to 
exhibit  di-morphism. 

Multicomponent  f luorozirconate  glasses,  a  relatively  new  class  of  materials, 
were  reported  first  by  M.  Chanthanasinh  (4)  in  1976  and  later  by  Poulain  et  al. 
in  1977  (5).  3ot.h  were  concerned  with  the  ZrF^-ThF^-Safi  system;  that  system 
in  also  the  subject  of  this  paper.  The  preparation  technique  typically  used 
by  earlier  workers  consisted  of  reacting  highly  purified  components  at  800  to 
900°C  cn  an  argon  atmosonere.  The  starting  materials  were  contained  in  a  her¬ 
metically  sealed  nickel  crucible  and  the  glass  was  formed  bv  pouring  the  full-.- 
reacted  melt  into  a  mole  residing  in  a  nitrogen-filled  glove  box.  The  vitreous 
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domain  for  ZrF—ThF^-Bar  can  be  described  as  an  area  or  a  triangle  sn  rbe 
ternary  diagram  bounded  by  the  following  points:  oo  mole  ”,  ZrF^,  ii  mole  1 
ThFq,  and  33  mole  Z  3aF-,.  RAP  preparations  rnroughouc  this  area  readily  yield 
transparent  glasses;  however,  the  composition  of  interesc  for  tnis  paper  is 
oOli  Ir?^,  71  ThF^ ,  and  33*  SaF-i,  which  is  essentially  the  center  of  me  vitre¬ 
ous  domain. 

This  paper  describes  the  RAP  approach  to  preparing  this  fluoride  glass 
utilicing  open  crucibles  and  controlled  atmospheres  of  anhydrous  HF  and  CC1.- . 
Physical  properties  including  spectral  data  are  given  as  well  as  glass  molding 
techniques . 


Experimental 


Purification  of  Materials 

ZrFj,  of  99.5.”  purify  was  obtained  from  the  Cerac  Company.  This  raw 
material  was  treated  with  dry  HF  at  400° C  and  sublimed  once  at  900° C  in  a 
stainless-steel  sublimator.  The  high-purity,  colorless,  transparent  needles 
were  then  combined  with  ThF.  and  3aF->  to  give  the  desired  ratio  —  0.60  Zr : 

0.07  Th:  0.33  3a.  The  ThF^  was  derived  from  99.9”  ThO-.  (Cerac).  The  oxide  was 
first  reacted  with  a  49”  HF  solution,  and  then  the  water  was  evaporated  away  at 
about  100oC  until  a  dry  powder  resulted.  :iext,  this  powder  was  slowly  heated 
and  subjected  to  gaseous  HF  uncll  melting  occurred.  The  high-purity  product 
was  the  source  of  Thr^  for  glass-forming  experiments.  The  Ear;,  from  EM 
laboratories ,  was  99.991  pure.  It  was  created  with  dry  HF  to  1000°C  prior  to 
use  as  a  starting  material  for  glass  preparations. 

Each  glass  forming  run  consisted  of  a  93-g  batch.  Vitreous  carbon  or  plat¬ 
inum  crucibles  were  used  in  an  Astro  Industries  (graphite  resistance;  furnace. 
The  furnace  was  vacuum  pumped  to  BOO^C  before  the  start  of  RAP.  The  reactive 
atmosphere  was  either  anhyarous  HF  diluted  to  10  mole  Z  with  helium,  or  gase¬ 
ous  CCI4  (Mailinkxodt  AR)  diiuced  with  helium.  The  mixture  was  neated  to 
900',C  in  5  hr;  the  resulting  melt  was  soaked  for  3  hr  and  chen  cooiea  abruptly 
to  room  temperature.  The  ingoc  could  then  be  annealed  at  ZZ0°C  to  remove 
strain. 


Class  Molding 

Molding  was  accomplished  with  a  10-ton  PHI  hydraulic  press  with  rsadaoie 
pressure  to  2.5  psl  intervals.  A  1/2-in.  cylindrical  tungsten  :arolde  die  was 
used.  The  inner  faces  of  the  plungers  were  polished  opcicaliy  fiat  co  i,  2 
fringe  (sodium  spectral  line  —  5396  A).  The  die  with  sample  inside  is  placet 
in  the  press.  The  cemperature  is  raised  to  HZ^C  and  increasing  oressure  is 
applied  until  a  pressure  drop  is  indicated  due  to  glass  flow.  Temperature  is 
held  tonstant  untii  the  die  is  filled  as  indicated  by  no  further  oressure  cron 
even  as  higher  pressure  is  applied.  The  temperature  is  then  lowered  below 
3G0°Z,  the  pressure  is  released,  and  me  sample  is  cooied  to  room  temperature. 

Results  and  Discussion 


Glass  Properties 

Early  batches  or  0.60  ZrF^tO. 37  ThF^:0.  33  3aF ->  glass  were  prepared  exclu¬ 
sively  with  a  RA  consisting  of  HF  diluted  with  helium.  Although  .'{-r.iv  powaer 
oactera  snowed  that  the  specimens  were  indeed  glassy,  they  ippeared  to  consist 
of  two  phases  —  a  primary  transparanc-coloriess  pnase  and  a  relatively  small 
cuancify  of  an  opaque-black  phase.  Electron  aicrcprobe  analysis  of  the  samples 
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indicated  chat  Che  opaque-olack  substance  contained  only  one-half  the  fluorine 
of  the  transparent  phase.  This  fluoride-deficient  phase  formed  regardless  of 
whether  a  platinum  or  carbon  crucible  was  used,  and,  surprisingly,  its  con¬ 
centration  appeared  to  increase  with  longer  HF  soak  times. 

The  formation  mechanism  of  thi3  opaque-black  phase  is  not  understood;  it 
would  appear  to  result  from  the  dissociation  of  ZrF4.~  Figure  1  shows  the  IF 
transmission  curves  (Beckman  IS  12)  of  three  fluoride  glass  samples.  Curve  A 
is  for  the  glass  0.60  ZrF4:0.10  ThF4:0.3Q  3aF->.  This  glass  is  manufactured  in 
France  by  La  Verre  Fluore  and,  using  an  inert  atmosphere,  is  prepared  with 
careful  attention  to  the  prevention  of  hydrolysis.  The  absorption  centered  at 
about  2.35  urn  may  be  the  result  of  a  trifling  amount  of  OH-  contamination. 

Curve  3  results  from  a  HF-prepared  glass,  0.60  ZrF4:0.07  ThF_, ;  0. 33  3aF->.  This 
specimen  appears  to  have  even  a  higher  concentration  of  CH~  than  does  A.  HF  is 
not  an  effective  reactive  atmosphere  for  these  glasses.  It  can  displace  OH- 
and  0“  from  simple  mixtures  of  3ar->  and  Thf4,  but,  at  temperatures  above  530°C, 
ZrO,  resists  conversion  to  ZrF4  (i.e.,  ZrOo  is  thermodynamically  favored 
over  ZrF4  even  in  the  presence  of  HF).  0H“  impurity  arises  from  hydrolysis 
stemming  from  water  emanating  from  the  interior  surfaces  of  the  furnace,  and 
HF  of  course  cannot  eliminate  this  water.  The  reduced  transmission  throughput 
is  caused  by  some  precipitated  fluoride  deficient  base.  Curve  C  is  0.60  ZrF^: 
0.07  ThF4:0.33  3aF2  glass  prepared  in  a  CCi4  atmosphere.  N’ote  that  the  trans¬ 
mission  is  continuous.  CCI4  is  an  extremely  effective  reactive  atmosphere.  It 
react3  with  both  OH-  in  the  melt  and  outgassed  water; 


CC1, 

*4 

CC1, 

>4 


(gas) 

(gas) 


20H , 


melt) 


+  2H,0, 

2  (gas) 


CO, 

CO, 


(gas) 

(gas) 


2HC1 ,  ,  +  2C1,  , 

(.gas)  (.melt) 


IHCd 


(gas) 


Moreover , 
chlorine, 
colorless 


CCi4  pyrolytically  cracks  at  temperatures  above  u00°C  to  give  nascent 
which  compensates  any  fluoride  deficiency  and  yields  remarkably 
and  transparent  Ingots. 


The  absorption  coefficients  of  sapphire,  Le  Verre  Fluore  glass,  and  Hughes 
f luorozirconate  glass  are  compared  as  a  function  of  wavelength  in  Figure  2. 

The  measurements,  made  with  a  Cary  Model  HR,  3how  the  Hughes  glass  1 CCl^  pre¬ 
pared)  to  be  less  absorbing  than  the  Le  Verre  Fluore,  but  not  nearly  as  good 
as  sapphire.  Fluorozirconate  glass  prepared  in  anhydrous  HF  gives  about  t.-.e 
same  result  as  Le  Verre  Fluore  glass.  The  physical  properties  of  the  glass 
0.60  ZrF4:0.07  ThF4:0.33  3aF2  prepared  an  CCi4  are  given  in  Table  1.  The  glass 
transition,  crystallization,  and  fusion  temperatures  were  determined  for  a 
sample  sealed  at  room  temperature  in  a  platinum  capsule  at  a  pressure  of 
<10-:3  Torr.  'later  insolubility  was  determined  by  soaking  typical  glass  samples 
in  tap  water  for  24  hr  at  room  temperature  and  ooserving  no  weight  chance. 


Molded  I lass 

Figure  3  3hovs  the  fluorozirconate  glass  (prepared  in  CC14)  and  the  die 
before  and  after  hot  pressing.  Pressing  was  done  at  312°C  and  1920  psi  in  air. 
Mote  che  frosty  areas  corresponding  to  flat  cut  corner  bevels  on  the  original 
blank.  Total  height  reduction  was  31.22  —  from  0.213  in.  to  0.150  in. 


Residual  gas  analysis  of  crystalline  ZrF^  shows  the  evolution  of  fluorine 
at  T  —  oOCC,  i.e., (white)  ZrF4  —  (black)  ZrF ,  —  F->. 
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WAVELENGTH  IN  MICRONS 

2.5  3.0  3.5  4.0  5.0  6.0  7.0  3.0 


WAVENUMBER,  CM"1 

FIG.  I 

Transmission  spectra  of  fluorozirconate  glasses.  Boch  faces 
polished  to  a  window  finish  for  all  samples. 


FREQUENCY  (WAVENUMBERS).  103  CM'1 
fig.  : 

UV  absorption  coefficient  of  uv  grade  sapphire,  La  Verre  Fluore 
and  Hughes  fluorozirconate  glass  vprepared  in  CCla  atmosphere). 
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TABLE  1 

Physical  Properties  of  0.60  ZrF^O.O?  ThF^:  0.33  Bax  2  Glass 


(Prepared  in  CC1,  Atmosphere) 

Properties 

— 

Values 

Optical  transparency 

0.3  um  GV  to  7  _m  IR3 

Density 

4.8  g/ca^ 

Glass  transition  temperature^ 

295®C 

Crystallization  cemperature^ 

390°C 

Fusion  temperature^ 

485  °C 

Refractive  index  (n^) 

1.53 

b  —2 

Knoop  hardness  (Kg-mm  ) 

250 

Humidity  effects 

Water  insoluble 

Heat  capacity  (C?)  at  45*0^ 
Coefficient  of  linear  expansion 

0.511  J/g-°C 

30  to  60°C 

i  * 

4.3  x  10  /“C 

250  to  270’C 

I 

. 

13.3  x  10'Vc 

1  ^ 

I  Rupture  strength 

( 

1 

j  35,300  psi 

1 

|  Absorption  coefficient  at  HF  laser 
wavelength  (2.8  urn) 

i 

s  .rt-3  -1 

b  x  10  ca 

i 

Absorption  coefficient  ac  DF  laser 

1  wavelength  (3.3  -m) 

1 

j  2  x  10  ^  cm 

1 

j  ^Prior  art  fluoride  glass  shows  a  st' 

:ong  absorption  at  2.35  _.n. 

Single  crystal  calcium  fluoride  (Car.,)  is  155. 

'  CPorzed  CaF.  is  14,000  psi. 

°Steasured  with  DuPont  thermal  analyser  model  1C90. 
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(a)  FACE  OF  DIG  PLUNGER 


(bl  FACE  OF  PRESSED  GLASS 


Hoc  pressed  Q.eC  ZrF  :0.0T  thr  i:';.  22  SaF  - 
glass  (prepared  »n  ICI^  acno30.--.ere  ■  .  2otr. 
surfaces  pnocograpre s  in  sod^ua  lig.ec 
while  sonracring  i  1.21  •  ootical 
reference  fl 


,ac  . 
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? igura  4(a)  snows  che  flatness  or  the  die  plunger  to  be  aoout  1/2  fringe  of 
sodium  Light  (5396  A/4).  Figure  4(b)  shows  chat  the  added  glass  surface 
closely  duplicates  the  die  surface. 


Siimiar' 

The  key  to  maximizing  IS  transparency  and  aechanical  strength  for  these 
halide  glasses  lies  in  eliminating  anionic  impurities  such  as  OH"  and  0“  which 
enter  the  condensed  phase  through  hydrolysis  during  che  mixing  and  melting 
stages  of  preparation.  The  HF  scrubbing  is  not  effective  in  this  case.  HF 
will  not  react  with  outgassed  water  emanating  from  che  furnace  interior,  anc, 
at  elevated  cemperatures,  it  will  not  react  with  ZrO->  to  yield  OrF^.  Moreover, 
the  use  of  HF  appears  to  promote  significant  quantities  of  a  black,  fluoride- 
deficient  phase.  CCl^,  on  che  other  hand,  eliminates  oxygen-containing  spe¬ 
cies  by  reacting  boch  with  atmospheric  water,  and  with  oxygen  species  in  che 
melt. 

Furthermore,  utilization  of  CCl^  for  scrubbing  compensates  fluoride  de¬ 
ficiency  by  adding  Cl~  to  the  F”  vacancies.  Fluoride  glass  prepared  with  the 
benefit  of  CCI4  RAP  is  continuously  transparent  from  0.3  urn  to  0.7  _m,  is 
unusually  hard  and  strong,  and  can  be  easily  molded  at  moderate  temperature  and 
pressure. 
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Optical  Properties  of  the  Surface  Piasmon  of  n-r  Silicon 

?.  F.  Robusto  and  R.  Brauns tein 

Department  of  Physics,  University  of  California,  Los  Angeles,  California  900 

The  optical  properties  of  the  surface  piasmon  of  n-F  silicon 
vere  determined  by  a  least  squares  curve  fit  of  the  Drude  intra¬ 
band  model  to  the  AIR  reflection  spectra.  The  surface  piasmon 
resonance  of  n+  silicon  is  highly  damped  (-u^T  -  2.69),  and  causes 
the  ATR  dispersion  cur/e  calculated  on  the  real  frequency  axis  to 
exhibit  backbending  for  its  complex  frequency- like  solutions', 
whereas,  the  theoretical  dispersion  curve  exhibits  the  standard 
vaveveccor  divergence  3t  the  surface  piasmon  frequency.  This  dis¬ 
crepancy  is  rectified  by  converting  the  optical  properties  of  n— 
silicon  into  functions  of  complex  frequency,  and  using  them  to 
calculate  the  ATR  dispersion  curve  in  the  complex  frequency  plane. 

A  comparison  of  the  Drude  dielectric  function  as  a  function  of 
real  and  complex  frequencies  was  analyzed.  The  comparison  shows 
that  the  difference  between  their  values  is  small  at  high  frequencies 
and  aeviate  the  largest  in  the  region  near  the  surface  piasmon 
resonance . 


Qptische  Messungen  des  Qberfldchen-Plasmons  von  n+Silizium* 

+ 

P.F.  Robusto  und  R.  Braunstein 

Physik-Abtei lung,  University  of  California,  Los  Angeles,  California  90024 

Die  optischen  Eigenschaften  des  Qberflachen-Plasmons  von  n+Silizium 
wurden  festgelegt,  indem  fur  das  Drude  Intraband-Model  1  die  Methode  der  kleln- 
sten  Quadrate  fur  das  Verfahren  der  stochasti schen  Kurvenermittlung  an  die 
abgeschwachten  Total reflexions-Spektra  Anwendung  fand.  Die  Qberflachen- 
Plasmonresonanz  von  n+SIlizium  1st  stark  gedampft  (^r  =  2,69),  und  bewirkt, 
daB  die  auf  der  wirklichen  Frequenzachse  berechnete,  abgeschwachte  Total  - 
reflexions-Dispersionskurve  Zuruckbiegung  fur  deren  komplexe  frequenzartige 
Losungen  aufweist,  wohingegen  die  theoretische  Di spersionskurve  die  Standard- 
Well  envektor-Divergenz  bei  der  Oberflachen-Plasmonfrequenz  aufweist.  Diese 
Unstimmigkeit  wird  korrigiert,  indem  die  optischen  Eigenschaften  von  n+  Sili- 
zium  in  Funktionen  komplexer  Frequenz  umgewandelt  werden,  und  indem  diese 
beniitzt  werden,  um  die  Kurve  der  abgeschwachten  Total  reflexion  in  der  Ebene 
der  komplexen  Frequenz  zu  berechnen.  Ein  Vergleich  der  dielektrischen  Funk- 
tion  nach  Drude  in  Abhangigkeit  von  wirklichen  und  komplexen  Frequenzen  wurde 
analysiert.  Der  Vergleich  ergab,  daB  der  Unterschied  zwischen  den  Werten  bei 
hohen  Frequenzen  niedrig  1st  und  daB  die  Werte  am  meisten  1m  Bereich  nahe  der 
Oberf lachen-Pl asmonresonanz  abweichen. 


I .  INTRODUCTION 


The  dispersion  relationship  of  surface  polaritor.s  cf  metals'1"’”  anc 
3 

semiconductors  has  been  extensively  studied  tor  materials  where  the  dancing 
of  the  surface  polariton  oscillation  is  small.  It  is  the  ain  of  this  study 
to  investigate  che  optical  properties  of  the  surface  piasmon  cf  n —  silicon 
for  which  the  oscillation  is  highly  damped.  The  general  dispersion  relationship 
for  surface  plasnons  for  an  air-material  interface  is"* 


ui  f  £  (to)  i 1/2 
c  ■+■  lj 


where  ki  i  is  the  surface  piasmon  wavevector  which  is  oarallel  to  the  inter- 

I  ■ 

face,  w  is  the  surface  piasmon  frequency,  c  is  the  speed  of  light,  and 
is  the  frequency  dependent  bulk  dielectric  function  of  silicon  which  for 
this  study  will  be  of  the  fora  of  the  Drude  free  electron  model  which  is 


o  (to)  =*  £,  (tu)  ■+■  i  sn(uj) 
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where  :  is  the  high  freauencv  dielectric  function,  o  =  •’ImNe” '>■;  m*)~  “  Ij 

CO  J  ^  to 

the  bulk  plasma  frequencv,  m*  is  the  electron  effective  mass,  is  the 
electron  concentration,  e  is  the  electronic  charge,  and  T  is  the  electronic 
relaxation  time. 

For  a  damned  surface  oscillation  the  ronrolex  frequencv  -  real  wavevector 
solutions  of  the  dispersion  relationship  given  by  Uq.  (1)  is  solved  by  me 
methods  of  Gammon  and  Palik.  The  wavevector  k  .  is  assumed  real,  while  tie 
frequency  of  the  excitation  is  complex  and  equal  to 
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where  oj,  is  the  surface  plasmon  frequency,  and  is  related  to  the  surrace 
plasmon  damping.  A  plot  of  jj^  versus  k|  ;  yields  the  surface  plastron  dispersion 
curve.  Briefly,  Gammon  and  Palik3  methods  involve  substituting  Eqs.  (1)  and  (3) 
into  Eq .  (1)  and  solving  for  a  -j,,  value  for  a  particular  -j,  value  such  that 
the  imaginary  part  of  Eq .  (1)  is  equal  to  zero.  With  these  values  ot  and 
n,  the  real  part  of  Eq .  (1)  is  solved  for  the  real  wavevector  k, i .  Repeating 
this  produces  a  set  of  a)1  ,  uj0  and  k|  i  values  are  obtained  from  which  the 
dispersion  curve  is  obtained.  The  complex  wavevector  -  real  frequency  solutions 
of  the  surface  plasmon  dispersion  relationship  [Eq.  (1)]  will  not  be  discussed 
in  this  study. 

The  Otto  ATR  technique1  is  a  method  for  coupling  incident  electromag¬ 
netic  waves  to  surface  polaritons.  3v  computer  curve  fitting  the  ATR  spectra 
to  the  Drude  dielectric  function  [Eq.  (1)]  the  optical  properties  of  the 
surface  active  medium  are  obtained.  Using  these  optical  constants  the  dis¬ 
persion  relation  of  the  surface  plasmon  on  the  real  frequency  axis  is  obtained 
from  the  computer  calculated  minima.  Each  minimum  is  obtained  by  fixing  the 
incident  angle  and  varying  the  frequency,  where  the  reflected  intensify  of 
p-poiarized  lighc  is  greater  than  90%  at  the  minimum. 

Section  II  contains  the  experimental  details  of  sample  preparation  and 
surface  plasmon  resonance  measurements.  Section  III  gives  the  results  of  the 
optical  measurements  and  discusses  the  dispersion  of  the  surface  plasmon  in 
the  real  frequency  axis.  In  Section  IV  we  discuss  the  ATR  anaivsis  m  the 
Comdex  frequency  piane  and  the  complex  frecuencv  optical  oroperties  ot  n-1- 
silicon.  In  Section  V  we  summarize  our  results. 
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ll.  experiment 

A.  Sample  Preparation 

The  n+  silicon  samples  were  obtained  by  diffusing  phosphorous  into  a 

5  ..-cm  n-cvpe  silicon  wafer.  The  initial  wafers  were  two  inches  in  diameter 

and  18  mils  thick,  and  were  polished  on  both  sides.  The  samples  were  orepared 

at  the  Micro-Electronics  Production  Laboratory  at  Hughes  Aircraft  Conpanv  in 

Carlsbad.  The  wafers  were  initially  prefurnaced  cleaned  in  a  50/50  mixture 

of  HnS0,  and  to  remove  anv  oxide  on  the  surface  and  to  take  off  anv  wax 

_  *1  _ 

left  from  the  polishing.  The  samples  were  then  diffused  at  925!5C  in  phosphine 
(?H^)  gas  for  10  minutes  with  a  preheat  and  post-heat  for  ten  minutes  in 
M,/0,.  After  the  n-1-  deposition,  the  samples  were  heated  at  1C00°C  in  pure 
nitrogen  for  130  minutes.  The  purpose  was  to  both  electrically  activate  the 
phosphorous  ions  and  to  diffuse  the  phosphorous  ions  into  the  bulk  of  the 
sample  ana  remove  the  high  concentration  of  ior.s  left  on  the  surface.  The 
samples  were  then  deglazed  in  a  10  to  1  solution  of  H.,0  and  HP  for  five 
minutes  for  the  purpose  of  removing  any  grown  oxide  on  the  silicon  surface. 

The  final  resistance  of  the  diffused  laver  was  o.l/scuare. 

The  initial  wafers  had  a  bow  of  0.115  mils  and  a  taper  of  0.0J  mils 
across  the  two  inch  diameter.  The  processed  wafers  were  scribec  and  or- ken 
into  ATS  sample  sizes  of  ifO  x  300  mils.  Aluminum  soaeers  were  electron  -earn 
evaporated  onto  two  comers  of  the  silicon  and  nine  '•  10  mil  diameter  ol uni rum 
spots  were  evaporated  between  the  two  corners.  The  silicon  wafers  were 
mounted  with  double  backed  tape  onto  optical  flats  12"  thick;.  \’:th  the 
silicon  d-'e  mounted  on  the  optical  flat  as  thev  were  ocmnressed  against  the 
base  of  the  prism  bowing  of  the  wafer  wculd  take  olace  when  onlv  the  two 
;orners  of  the  lie  were  coated  with  aluminum.  lonsepuentl" .  the  audit:  -r. 
the  ••  20  mil  diameter  spots  across  the  center  region  reduced  t:.<‘  wine 
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signif i candy  and  allowed  a  reasonably  uniform  air  gap  to  exist  between  the 
prism  base  and  the  silicon. 

3.  Experimental  Measurements  and  Data  Analysis 

Initial  investigations  of  the  optical  properties  of  rm  silicon  were  done 
by  measuring  the  near  normal  incidence  reflectivity  of  silicon  versus  gold. 

These  measurements  were  not  absolute  but  were  done  for  the  purpose  of  finding 
the  plasma  edge  and  minimum.  The  minimum  appeared  to  be  at  1.5  _  and  conse¬ 
quently  ATR  spectrum  runs  were  performed  in  che  range  of  2  t  to  9  - .  AT?, 
measurements  were  performed  by  measuring  the  ratio  of  the  reflected  intensity 
of  p-poiarized  light  (R  )  with  the  sample  pressed  against  the  ATR  prism  to 
that  with  the  sample  removed.  The  Otto  ATR  method  of  exciting  surface  plasmon 
modes  was  used  where  the  surface  plasmon  was  for  an  air-Si  interface. 

A  Caf,  prism  (37. 1”  -  65. 6’ [apex]  -  57.2°)  was  used  for  all  measurements; 
it  has  an  average  index  of  1.45  and  it  transmits  out  to  9  microns.  The  prism 
base  was  flat  to  greater  than  a  tenth  of  a  wave  (0.5461  u)  and  the  uniformity 
of  the  air  gap  was  determined  by  viewing  Newton  interference  fringes  for 
ingles  less  than  the  critical  angle  and  by  viewing  the  surface  plasmon  resonance 
of  aluminum  which  appears  in  the  visible  spectrum  for  angles  greater  than  the 
critical  angle,  uniformity  of  this  color  absorption  was  adjusted.  7cod  uni¬ 
formity  of  the  air  gao  across  the  whole  silicon  surface  was  impossible  to  obtain 
hue  to  its  natural  bow;  consequent!’/ ,  the  innut  beam  was  masked  to  hit  :t’.T;  the 
uppr oximateiy  uniform  section.  2ue  to  this  bowing  orcoiem,  uniform  air 
:ouiu  be  obtained  only  for  small  spacings  between  the  silicon  sample  and 
the  prism  base.  The  probable  reason  for  this  reduction  of  bowing  at  small 
spacings  is  that  both  nigh  points  on  the  silicon  and  dust  ■'.articles  are 
present  between  the  prism  base  and  silicon,  surface  which  acts  as  a  no: ura. 
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The  globar  light  source  was  focused  down  f/oeam  civergence  -  5°)  and 
optical  alignment  of  the  beam  could  be  seen  by  viewing  the  output  beam  of 
the  prism  assembly  with  a  mirror.  The  reflection  spectra  were  taken  on  a 
Perkin  Elmer  model  301  far  infrared  spectrometer  operated  in  a  single  beam 
mode  with  the  addition  of  appropriate  gratings  and  filters.  The  13  cps  chopped 
signal  from  the  thermocouple  detector  was  digitized  and  then  computer  av-ruged. 
The  wire  grid  polarizer  was  placed  in  che  p-ooiarizing  position  and  placed 
before  the  input  beam  to  the  prism. 


C .  Optical  Constants  and  ATR  Curve  Pitting 

The  ATR  reflection  spectra  were  taken  at  room  temperature  and  displayed 

in  Figs.  1  and  2.  The  spectra  are  plotted  as  normalized  reflectivity  ratio 

R  versus  wavelength  where  the  data  points  are  corrected  for  the  effects  of 
o 

the  aluminum  spacers.  The  optical  properties  of  the  r.-r  silicon  were  determined 
by  doing  a  least  squares  curve  fit  of  the  crude  parameters  M/m*  anc  and  the 


air  C30  d.  The  high  frequencv  dielectric  constant  c  used  in  the  Truce  model 
Eq .  (2)  was  taken  to  be  equal  to  ^  =»  11.7  the  dielectric  constant  of  silicon 
at  zero  frequency  in  the  absence  of  free  carriers.  iince  silicon  is  honepoi^r , 
the  iieiectric  constant  or'  silicon  is  equal  to  c  in  the  absence  of  tree 


carriers  from  zero  frequency  up  to  the  first  interbund  transition.  The  iteiectr: 


instant  c 


Li.  7  was  determined  by  Suizberg  anc  -’ilia  on  measuring  the 


refractive  index  of  a  prism  made  of  single  crystal  silicon.  This  dielectric 
constant  also  agrees  with  the  value  measured  at  raui;  frequencies  :n  intrinsic 
silicon  by  Dunlap  and  Valters.' 

The  parameters  M/m*  and  -  were  varied  sior.g  with  t ne  air  gar?  a  to  obtain 
the  best  fit  least  squares  parameter  aefinec  as 


.  ;<r  -  <?.  r 

1  - ->  inner;;  o  ext 
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RESULTS  AND  DISCUSSION 


A.  Optical  Properties 

The  Drude  free  electron  model  of  che  dielectric  function  fits  sacis- 
factorily  to  the  ATR  soaccra  curves  as  shown  in  Tigs.  1  and  2.  The  experi¬ 
mental  error  is  relatively  large  due  to  the  bowing  of  the  silicon  chip  and 
the  need  to  add  the  extra  aluminum  spacer  pads.  The  deviation  in  the  region 
of  the  minimum  for  the  curve  fit  in  Tig.  2  is  possibly  indicative  the  need 

to  include  a  frequency  dependent  collision  time.  A  similar  deviation  was 

9 

obtained  by  Kukharskii  and  Suboshiev  when  they  tune  fitted  a  Drude  model 
to  n+  silicon.  They  corrected  this  deviation  by  using  a  frequency  dependent 
relaxation  time.  There  exists  other  problems  in  this  part  of  the  spectral 
region  where  the  Drude  model  should  be  corrected.  N'aneiy,  as  shown  by  Saiobe 
ana  Villa,  increases  from  11,7  at  p  c  to  '•  11.32  at  2.3  u.  This  Is  duo 
to  the  dispersion  of  the  index  of  refraction  when  the  band  edge  is  approached 
opitcer  and  Jan*^  have  shown  that  this  is  a  small  absorption  band  in  the 
range  of  1.5  to  5  d  which  is  dependent  on  carrier  concentration.  Ccrrectic 
for  these  effects  were  not  included  due  to  the  relatively  large  experimental 
error  and  in  the  region  '<  >  5  -  the  region  where  the  surface  daemon  exist  • 
the  curve  tits  were  good  enough  and  should  accurately  reproduce  the  surraoi.- 
nlasmon  response. 


The  experimental!'/  iecemir.ed  Drude  free  warrior  -artmeters  a 

faole  I  along  with  the  opcicai  cicbilcry  .  =  e- 'm*  and  the  electron 

xcncentrat ion  ".  The  overall  sc curacy  or  this  least  squares  fitti 

an  be  seen  bv  viewing  Tig.  2  where  each  parameter  :  s  varied  '-v  12 

values,  f.  -m  the  best  fit  which  is  also  clotted  wi tn  it.  A  value 

,ya  3 

N  ’  '  •  13*'  e.  :m  was  ietermcr.ee  from  a  measurc-m«?r.t  : :  electrics. 
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(.1.32  <  10  2-cn)  and  using  the  resistivity  versus  impuritv  ::ncer.crat:on 

?  1 

listed  in  Sze--:  this  agreement  is  reasonable  and  the  difference  is  probablv 

due  to  ditterenc  sample  mobilities.  Kukharskii  anc  Suboshiev'  obtuinea  a 

-15 

relaxation  tine  of  11.27  <  10  sec  tor  an  n-cype  silicon  or  carrier 

19  3 

concentration  5.5  *  10  '  e/cm  .  Thev  also  obtained  a  value  or  a  T  =  2.69 

p 

which  is  in  agreement  with  the  present  results. 


3.  Disoersion  Ar.alvsis 


The  value  of 


.69  as  determined  from  the  Orude  ootical  naraneters 


(Table  I)  shows  that  the  surface  plasmon  resonance  of  rH-  silicon  is  fairly 
highly  damped.  The  complex  frequency  -  real  wavevector  solutions  of  the 
theoretical  dispersion  curve  Eq .  (1)  is  shown  in  Figs.  1  ana  3.  The  complex 
frequency  solutions  exhibit  the  standard  behavior  of  the  surface  plasncn 
resonance;  namely,  and  a.,  •*  -  1/2"  as  the  wavevector  k  -  «. 

A.  5  - 

The  AT?,  surface  plasmon  dispersion  curve  is  found  by  plotting  the  energy 
and  wavevector  at  the  minimum  reflectivity  point  in  the  AT?  spectrum.  The 
AT?  minimum  point  is  taken  when  the  reflectivity  of  p-poiariied  light  is  at 
90.1 .  The  general  computer  procedure  has  been  previouslv  explained,*-  and  the 
results  of  the  analvsis  are  shown  in  Table  17.  Ir.  this  table,  vou  tan  see 


that  as  the  angle  AO  angle  -wit1’,  reset 


tne  prism  case'  is  increases, 


wavelength  at  the  minimum  '.L  increases  to  longer  wavelengths  and  smaller 


exoer  i  meat  a  -  wavevec  cor ; 


.are  acn levee.  . .us  ar.a.v? 


the  incident  angle  is  decreased,  the  AT?,  listers icn  curve  ;p-r caches 


-me,  as  t.i r.  :e  seer 


ituavir.g  the  last  column  wnior.  gives  the  '.at;  ; 


wavevector  generate-!  •••  tne  T'.sa 


to  the  wavevector 


..il.  7:-.. it  1  i ,  the  ettrerimc-r.t.al  rati- 


■a it n  r  tac ::<• 
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lighc  line  and  then  back-bends  going  into  3rewster  nodes  where  k  ■  ■'  c . " 
Therefore,  the  ATS  determined  dispersion  shows  no  resonance  as  is  achieved 
bv  the  theoretical  dispersion  relationship.  Since  c  ~  =1.69,  there  is  no 

p 

physical  explanation  why  there  is  no  resonance  present  in  the  ATS  dispersion 
curve  obtained  using  real  frequencies.  Consequently  the  only  true  coupling 
to  the  surface  plasmcn  resonance  of  n+  silicon  will  occur  by  probing  the 
complex  k  response  which  is  experimentally  achieved  by  fixing  the  excitation 
frequency  and  by  varying  the  angle  of  incidence.  The  complex  •*>  response 
measurements  will  only  allow  moderate  coupling  which  is  achieved  by  riding 
on  the  ridges  and  valleys  of  the  surface  plasmon  excitation  in  the  real 
frequency  plane. 
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17.  ATR  ANALYSIS  IN  THE  COMPLEX  FREQUENCV  FLA. NT  .AND 
COMPLEX  FREQUENCY  OPTICAL  PROPERTIES  OF  IH-  SILICON 

Due  co  che  discrepancy  between  the  theoretical  damped  surface  plasmcn 
dispersion  relationship  which  shows  a  definite  resonance,  and  the  ATR  minima 
dispersion  curve  which  shows  back-bending  which  is  calculated  using  real 
rrequency  dependent  optical  constants,  an  aporoach  was  made  to  calculate  the 
ATR  reflection  spectrum  in  the  complex  frequency  plane.  In  order  to  calcu¬ 
late  che  ATR  spectrum  in  the  complex  frequency  plane,  che  optical  constants 
of  the  surface  active  medium  and  the  Fresnel  reflectance  coefficient  must  be 
redefined  to  be  functions  of  complex  frequency.  Since  the  frequency  is 
complex,  then  the  real  [-;.(_;)]  and  imaginary  ;  parts  of  the  dielectric 

functions  are  themselves  compiex  as  given  by 

e '•  3  C.  '.-)  ~  i  SV.u) 

where 

=  Re  (j)  -r  i  In  (u) 

-  1  X 

=  Re  +  i  In 

where  j  is  the  complex  frequency  defined  by  —  iw,.  If  one  substitute 


Eas , 

Cdb;  and 

(. ic'1  into  Eq.  1 -*a) 

and  regroup  int 

o  real  ana 

iT.a2ir.arn  -aris. 

one 

obtains 

:  a ) 

=  [Re  Z,(~)  -  Im  ; 

,'~)3  -  i  [Re 

~  In 

■1 .  « )  ’ 

- 

*■ 

*■ 

New 

using  the 

standard  definition 

of  £,  !  *.)  and  2 

% j  be 

equal  : e  rr.e  re  a 

and 

imaginary 

parts  of  the  dielec 

trie  function  r 

esrecr ivel 

*  as  :iven  :  r  be  . 

v  5 ) 

we  obtain 

fa)  =  i _ > , ,  -  Re  ,  „>  -  In  ^ j 

f~>)  -  »  Re  •  _j  -  Im 


Bv  redefining  -he  dielectric  fnncricr.  in  this  wav,  ye  show  that  the 
efface  of  the  coco  lex  frequency  is  to  nix  the  old  relationship  or  z.  and 
That  is,  now  che  real  part  of  the  dielectric  function  is  affected  by  -- , 
through  its  imaginary  rare.  Also  by  the  definition  of  the  retro  lex  frequency 
cnac  as  -*■  0,  the  imaginary  parts  of  C,  (-0  and  i0(u)  both  gc  to  tero  and 
we  have  the  original  relationship  in  the  real  frequency  plane.  ".v'ith  the 
dieieecric  function  as  defined  in  Bqs.  (03}  and  'Voi,  we  can  use  the  standard 
relationships  between  the  jptical  conscancs  n  and  k  and  the  dielectric  func¬ 
tion  £,(-)  and  £,(-)  in  the  complex  frequency  plane  which  is 


n  =■  =  i,  4  [  (i,  " 


’) 


rC  -  , -J-J 


Consequently  the  complex  index  on  refraction  tf  the  medii 
unction  of  :o npiex  frequency  and  mar  be  used  directly  in  net 
'eflectivitv  in  the  complex  frequency  plane.  This  is  dene 
‘resne.  reflection  coefficient  for  o-noi ar i cec  licht  as: 


:  e  a  e : i n i n £  one 


cos%  -  cos?. 

- _ 0 _ j _ - 

cosd  ^  ^  cos? 


:os?.  - 


•cos-. 


?  3  T 


where  h'„  is  the  complex  index  of  refraction  of  the  surcace  active  me 

j 

in  this  case  is  silicon,  and  ",  is  the  ccmoiax  inaex  0:  refraction  : 
material  between  the  prism  and  the  surface  active  meuium.  wnicr.  is  0 
case  study,  and  7 ,  and  i,  are  the  incident  angles  for  meuium  1  ar.c  . 

—  w 

The  conn  lax  frequency  deper.den  e  jr  the  real  md  inagir.arv  •'art 
index  of  refraction  can  also  be  -bt aimed  using  a  iif  terer.c  -'r 
method  is  more  difficult  tv  solve  and  involves  substituting  iq  s • 
into  Bos.  1  ~ \i  ard  ~b.;  and  final."  r ’defining  the  r.-mr :  no-.-  '■ 


is 
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N  =  n  +  ik 

where 

n  =•  nCo^.ci^)  =  Re  n  -  Im  k 
k  *  kCuj^.w^)  “  Re  k  +  Im  n. 


(9a) 


(9b) 

(9c) 


Implementing  this  procedure,  the  optical  constants  of  n+  silicon  were 
expressed  as  functions  of  complex  frequency  and  used  to  calculate  the  ATR  dis¬ 
persion  curve.  The  calculations  were  done  by  using  a  high  speed  computer  pro- 

g 

gram.  The  general  procedure  was  to  first  obtain  a  set  of  anc*  ^  values  as 

determined  from  calculating  the  complex  frequency  -  real  wavevector  dispersion 

curve  [Eq.  (1) ] ;  that  is,  by  solving  the  imaginary  part  of  the  dispersion 

relationship  a  value  for  jj  was  obtained  for  each  u,  value.  This  set  of  jj,  and 

l  L  1 

values  was  used  to  calculate  the  real  and  imaginary  parts  of  anc*  £->  from 
Eqs.  (2)  and  (4);  and  then  using  these  values  in  Eqs.  (6a)  and  (6b)  to  obtain 
Sl^i)l,ai2^  an^  C2^W1,W2^’  are  used  in  Eqs.  (7a)  and  (7b)  to  obtain 

the  complex  frequency  dependent  optical  constants.  These  complex  frequency 
optical  constants  and  Fresnel  reflection  coefficients  were  used  to  calculate 
the  ATR  dispersion  curve  in  the  same  way  as  done  in  the  real  frequency  spectrum; 
but  now  the  ATR  minimum  is  found  by  fixing  the  angle  of  incidence  and  by 
scanning  through  uj^  the  real  part  of  the  surface  plasmon  frequency. 

The  ATR  dispersion  curve  minima  are  plotted  in  Fig.  6  along  with  the  damped 
surface  plasmon  dispersion  curve.  There  is  no  discrepancy  between  the  ATR 
dispersion  curve  calculated  in  complex  frequency  plane  and  the  theoretical 
dispersion  curve.  The  complex  frequency  ATR  analysis  shows  a  definite 
resonance.  The  error  bars  indicate  the  range  of  uncertainty  that  the  ATR 
minima  were  obtained.  The  deviation  at  lower  values  of  wavevector  is  due  to 
the  way  that  the  ATR  dispersion  curve  was  generated,  that  is,  with  fixing  the 
angle  of  incidence  and  varying  the  frequency  of  the  input  light.  This  type  of 
ATR  analysis  tends  to  ride  along  valleys  or  ridges  of  the  dispersion  at  lower 
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waveveccors. 


Using  the  complex  frequency  dependent  dielectric  function  the  surface 
plasmon  obeys  the  following  resonance  condition: 


where 


*  -  1 


£l('i°l’‘V  *  Re  £1  ~  Im  £2  3  "  l 

£.(w.,aj7)  »  Re  e,  +  In  ■  0  .  (101 

The  real  and  imaginary  parts  of  the  complex  frequency  dependent  Drude  dielec¬ 
tric  function  is  shown  in  Figs-  7  and  8.  As  can  be  readily  seen,  the  deviation 
of  the  real  and  imaginary  parts  of  the  dielectric  function  becomes  large  near 
the  surface  plasmon  resonance  point.  In  Table  III  there  is  a  comparison  of 
the  values  of  the  real  part  of  £ (oo)  which  is  a  function  of  complex  frequency 
and  £^(10)  calculated  from  Drude  model  for  real  frequencies.  As  can  be  seen 
from  Table  III,  the  difference  in  the  region  close  to  the  surface  plasmon 
resonance  (1  ■  6.66  ul  is  quite  large.  The  wavevector  solutions  of  the  surface 
plasmon  dispersion  relationship  [Eq.  (1)3  in  the  region  between  £(w)  *  -1  and 
e(ui)  »  0  are  imaginary  and  are  not  surface  modes  in  the  fact  that  they  are 
uncoupled  modes  between  the  incident  light  and  the  surface  active  medium. 

The  complex  frequency  dielectric  functions  shown  in  Figs.  7  and  8  have 
extrapolated  lines  covering  this  region.  Also  in  the  region  close  to  c(jj)  =  -1 
the  complex  solutions  are  very  difficult  to  obtain  due  to  the  nature  of  the 
resonance,  that  is,  the  solutions  occur  in  a  very  narrow  valley  with  a  step 
fall-off  on  both  sides  of  the  valley  and  can  be  easily  missed  by  the  computer 

Q 

procedure.  The  surface  plasmon  resonance  point  was  determined  and  is  marked 
by  an  X  on  the  dielectric  function  plots. 

In  Table  IV  a  comparison  is  made  of  the  dielectric  functions  obtained 
for  real  and  complex  frequency.  The  dielectric  functions  approach  one  another 
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ac  high  frequencies  and  deviate  the  largest  in  the  area  of  the  surface 

plasmon  resonance.  Solutions  are  shown  in  the  regions  where  e(o>)  >  0  and 

e(oj)  <  -1;  where  for  e Coj)  >  0  the  dispersion  analysis  yields  solutions  for 

Brewster^  modes,  and  for  s(ui)  <  -1  solutions  are  surface  polariton  modes 

(Fano  modes) .  The  surface  plasmon  resonance  occurs  when  *  jjg  and 

^2  m  -1/2X  where  T  is  the  electron  relaxation  time.  Solutions  for  the 
13 

Brewster  modes  show  that  the  bulk  plasmon  frequency  is  also  complex  and 

occurs  at  jJ,  »  ui  and  to.  =»  -  l/2x. 

1  P  2 


V.  CONCLUSIONS 


It  can  be  concluded  from  analyzing  the  surface  plasmon  of  n+  silicon 
that  the  dispersion  relationship  obtained  from  plotting  the  ATR  minima  obtained 
in  the  real  frequency  plane  does  not  exhibit  the  surface  plasma  resonance; 
whereas  the  theoretical  dispersion  relationship  yields  the  resonance. 

Consequently,  for  this  n+-  silicon  material  which  has  a  relatively  high 
damping  the  surface  plasmon  will  not  exhibit  a  resonance  in  the  real  frequency 
spectrum  when  coupling  is  performed  using  the  ATR  experimental  technique  of 
fixing  incident  angle  and  varying  frequency.  By  using  the  complex  frequency 
analysis  either  through  the  theoretical  dispersion  relationship  or  using 
our  developed  complex  frequency  dependent  optical  constants  and  analyzing  the 
ATR  minima  in  the  complex  frequency  plane  the  location  of  the  complex 
frequency  resonance  can  be  obtained. 

A  comparison  of  the  real  and  imaginary  parts  of  the  Drude  dielectric 
function  for  n+  silicon  for  both  real  and  complex  frequencies  shows  that 
they  numerically  deviate  the  largest  in  the  area  of  the  surface  plasmon  and 
approach  one  another  in  value  at  high  frequencies. 
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Figure  Captions 

Fig.  1.  Normalised  ATR  reflection  spectrum  of  n+  silicon  plotted  as 

normalised  p-polarized  reflectivity  versus  wavelength.  Solid 
points  are  the  data  points  corrected  for  effects  of  the  aluminum 
spacers.  Solid  curve  is  the  best  fit  least  squares  curve  with 
the  best  fit  parameters  shown  in  the  blocked  region.  Shown  in 
the  insert  is  the  ATR  geometry  used  with  the  angle  of  incidence 
of  che  input  beam  is  5  degrees. 

Fig.  2.  Normalised  ATR  reflection  spectrum  of  n+  silicon  plotted  as 

normalised  p-polarized  reflectivity  versus  wavelength.  Solid 
points  are  the  data  points  corrected  for  effects  of  the  aluminum 
spacers.  Solid  curve  is  che  best  fit  least  squares  curve  with 
the  best  fit  parameters  shown  in  the  blocked  region.  Shown  in 
the  insert  is  the  ATR  geometry  used  with  the  angle  of  incidence 
of  the  input  beam  is  0  degrees. 

Fig.  3.  ATR  reflectivity  spectrum  of  n+  silicon  where  each  of  che  best 
fit  parameters  T,  N/m*,  and  d  are  changed  to  15%  higher  values. 

Fig.  A.  Theoretical  dispersion  curve  [Eq.  (1)]  of  the  surface  plasmon  of 
n+  silicon  plotted  as  the  real  part  of  the  complex  energy 
versus  real  wavevector  k, j .  Also  shown  is  the  light  line  (k  *  jj/c) 
along  with  the  energy  wavevector  light  line  (dashed  line)  generated 
by  che  ATR  configuration  shown  in  Fig.  2.  The  theoretical  disper¬ 
sion  curve  is  plotted  from  Che  parameters  given  in  Table  I. 
Imaginary  part  of  the  complex  energy  versus  real  wavevector 

for  the  surface  plasmon  of  n+  silicon. 


Fig.  5. 
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Fig.  6.  Same  diagram  as  in  Fig.  4  with  the  inclusion  of  the  ATR  minima 

points  calculated  using  the  complex  frequency  dependent  dielectric 
functions  as  defined  in  Eq.  (4).  Error  bars  show  the  range  within 
which  the  minima  were  determined. 

Fig.  7.  Real  and  imaginary  parts  of  complex  frequency  dependent  versus 
wavelength  for  n+  silicon;  where  e,  (uj,  ,  lO  *  Se  e,  +  i  Im  s,  and 

■L  1  1  1  X 

is  not  redefined  as  in  Eq.  (4)  and  where  L.,  =  2Tc/u^. 

Fig.  8.  Real  and  imaginary  parts  of  complex  frequency  dependent  e7  versus 
wavelengths  for  n+  silicon;  where  =>  Re  e,  +  i  la  £,  and 

is  not  redefined  as  in  Eq.  (4)  and  where  *  2"Xc/bi^ . 


TABLE  I.  Experimentally  Determined  Drude  Free  Carrier  Parameters  for  n+ 


Silicon 


.(1) 


N/m*  T  N  u 

'  opt. 

(x  1020  e/cm3)  (x  10"15  sec)  m*^2)  (x  1019  e/cm3)  Ccm1 2/V- 


11.7 


3.58 


8.62 


0.26 


9.31 


58, 


(1)  Input  parameter  for  calculation 

(2)  Fraction  of  electron  rest  mass  obtained  from  Spitzer  and  Fan.^ 


sec) 


3 


•ABLE  II. 


m 


i'i) 

Analysis  of  ATR  Calculated  Dispersion  Curve  for  n+  Silicon' 


(2) 

A  0 

(degrees) 

L  (3) 

0 

(microns) 

D  C3) 

0 

(microns) 

■.exp 

<  i  i 

i  ] 

(X  104 

(4) 

cm  2) 

,exp,  (5) 
k||  /ko 

23.8 

468.9 

400.9 

1.35  x 

-0 

10  ~ 

1.00886 

17.4 

1074.4 

11.4 

1.08  x 

10“2 

1.839 

49.2 

1164.6 

10.97 

1.02  x 

10“2 

1.894 

53.1 

1376.7 

10.34 

9.126  x 

10~3 

1.9997 

57.1 

1685.4 

10.11 

7.325  x 

io"3 

2.099 

61.08 

2117.5 

10.24 

6.493  x 

io“3 

2.188 

66.88 

3195. 7 

11.182 

4.519  x 

10-3 

2.279 

(1)  Input  Drude  parameters: 

=  11-7. 

T  =  5.62 

x  10~15 

sec,  N’/m*  = 

20 

3.58  x  10  . 


(2)  Angle  of  incidence  with  respect  to  the  prism  base  where  n  =2.5. 

? 

(3)  and  are  the  wavelength  at  the  minimum  and  Che  air  gap  respectivitv 
to  obtain  a  p-polarized  reflectivity  of  90%. 

(1)  This  is  the  experimental  wavevector  at  the  minimum  where  k.^  = 

i  i 

P-  n  sin(AO) . 

Lo  ? 

(5)  This  is  che  ratio  of  the  experimental  wavevector  to  the  wavevector  of 


light  (k.  =  vaj/c  ■  2t/L  )  . 

o  o 


vjJ/  c 


I 


TABLE  III.  Comparison  of  Che  Real  Parc  of  the  Complex  u>  Dieleccric  Functions 


£,  (to)  and 

i 

£,((u)  versus 

the 

Dielectric 

Function  for 

Real  Frequency 

t 

jJ 

L(p)( 

1}  Re  sl(2) 

e1Cia') 

(3) 

-X<3> 

Re  £2C4) 

e2(uj') 

c 

wC.  n 

2 

10.4351 

10.4351 

0 

0.155894 

0.155894 

0 

3 

8.90629 

8.90623 

0.00006 

0.516462 

0.516484 

0.00002 

4 

6.85904 

6.85776 

0.001 

1.19291 

1.19358 

0.0007 

5 

4.39286 

4.37016 

0.023 

2.2433 

2.25844 

0.015 

5.5 

3.09752 

2.98949 

0.108 

2.87194 

2.95223 

0.080 

5.9 

2.31265 

1.82975 

0.419 

3.20381 

3.5886 

0.385 

6.66 

0.229251 

-0.486703 

0.716 

4.34137 

5.00156 

0.660 

6.7 

0.0409945 

-0.612146 

0.653 

4.4751 

5.08339 

0.608 

6.3 

-0.43728 

-0.927106 

0.4898 

4.8232 

5.29124 

0.^68 

7 

-1.28267 

-1.56256 

0.2799 

5.44094 

5.72098 

0.280 

7.5 

-3.06521 

-3.1795 

0.114 

6 . 75068 

6.87692 

0.126 

8.5 

-6 . 46159 

-6.50066 

0.039 

9.48172 

9.53347 

0.052 

9 

-8.15995 

-8.18639 

0.026 

10.990 

11.0292 

0.039 

10 

-11.5531 

-11.5672 

0.014 

14.3139 

14.338 

0.024 

11 

-16 . 5556 

-16.5623 

0.007 

20.0136 

20.0236 

0.015 

(1) 

Wavelength  of  t 

he  incident  light 

in  microns 

and  for  the 

complex  frequenc 

ies 

is  L  =■  2tc/uj.  . 

(2)  Real  part  of  che  complex  dieleccric  function  €^(uj)  =  +  i  jO  = 

uO  =  Re  +  i  Im 

(3)  As, 2  3 4  jRe  ,  jj9)  -  e^(w')1  where  and  are  the  real  and  imaginary 

parts  of  the  complex  frequency  jj,  and  cj'  is  the  real  frequency. 

(4)  Imaginary  part  of  the  complex  dielectric  function  r„(u)  »  ; 

=  Re  s,  *  i  Im  £,,. 

l£n  =  ']Re  uj^)  -  £„(«’)!  . 


(5) 


TA3LE  ;v.  Comparison  of  the  Drude  Dielectric  Function  as  a  Function  of  Seal 
and  Complex  Frequencies 


L(u) 

(1) 

£,(*’)< C2) 

m./3) 

£„(u>')(4) 

L. 

£2Coj)(:> 

2 

10.4351 

10.4343 

0.155894 

0.154233 

3 

8.9062 

3.90199 

0.516484 

0.501147 

4 

6.85776 

6.32434 

1. 19358 

1.10091 

5 

4.37016 

4.16524 

2.25844 

1.76509 

5.5 

2.98949 

2.50169 

2.95223 

1.73513 

5.9 

1.82975 

0.91369 

3.5686 

0.67388 

6 . 66 

-0.486703 

-1.91911 

5.00156 
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j  ^  +  ia.,  as  defined  in  £q.  (4),  that  is,  £.  fa)  =  s,  (a,  ,  a-,)  • 

(4)  Imaginary  part  of  dielectric  function  as  a  function  of  real  frequency  j' 

(5)  Imaginary  part  of  dielectric  function  as  a  function  of  complex  frequency 

-J  =  (*>  +  ia)0  as  defined  in  Eq .  (4),  chat  is,  £„(a)  =  a„)  . 
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Optical  Measurements  of  the  Surface  Plasmon  of  Copper 

?.  F.  Robusto  and  R.  Sraunstein 

Department  of  Physics,  University  of  California,  Los  Angeles,  California  9002^- 


The  optical  properties  of  the  noble  metal  copper  were 
determined  by  compucer  fitting  the  experimentally  obtained  attenu¬ 
ated  total  internal  reflection  (ATR)  spectra  of  its  surface  plasmon 
in  the  visible  spectral  region  (0.42  a  to  0.62  d) .  The  dispersion 
curves  of  the  surface  plasmon  resonance  of  copper  were  determined; 
they  show  that  the  resonance  exists  on  the  steep  absorption 
edge  of  S,  at  2.16  eV.  This  absorption  edge,  which  is  due  to  the 
d-oand  to  Fermi  level  transition,  causes  the  large  shift  in  the  loca¬ 
tion  of  the  surface  plasmon  excitation  from  i  U  )  *  -1  to  £,<w  }  * 

13  1  S 

-4.7  and  causes  the  true  resonance  to  occur  in  the  complex  frequencv 
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Optische  Messungen  des  Qberflachen-Plasmons  von  Kupfer* 

4. 

t  ~ 

P.F.  Robust o  und  R.  Braunstein 

Physik-Abteilung,  University  of  California,  Los  Angeles,  California  90024 

Die  optischen  Eigenschaften  des  Edelmetalls  Kupfer  wurden  festgelegt, 
indem  die  experimental!  erhaltenen,  abgeschwachten  inneren  Totalreflexlons- 
Spektra  (ATR,  abgeschwachte  Total  reflexion)  dessen  QberfVachen-Plasmons  durch 
Computer  in  den  sichtbaren  Spektralbereich  angepaQt  (0,42  p  bis  0,63  p)  wur¬ 
den.  Die  Dispersionskurven  der  Oberflachen-Plasmonresonanz  von  Kupfer  wurden 
bestimmt;  sie  zeigen,  daG  an  der  steilen  Absorptionskante  von  Resonanz  bei 
2,16  eV  vorhanden  ist.  Oiese  Absorptionskante,  die  auf  den  Obergang  vom 
d-Band  zum  Fertni-Niveau  zuruckzufuhren  ist,  bewirkt  die  groSe  Verschiebung  in 
der  Lage  der  Qberflachen-Plasmonanregung  von  (ujs)  =  -1  auf  (aj$)  =  -4,7 
und  bewirkt,  daG  die  echte  Resonanz  in  der  komplexen  Frequenzebene  auftritt. 


I.  INTRODUCTION 


The  optical  constants  of  copper  have  been  measured  since  the  time  of 

Drude.  Copper  is  known  to  exhibit  interband  transitions  throughout  the  visible 

region  of  the  spectrum  starting  with  a  d-band  to  Fermi  level  transition  at 

v  2.1  eV.  However,  copper  is  known  not  to  exhibit  a  bulk  plasma  resonance. 

The  general  explanation  of  this  is  that  due  to  the  d-states  the  positive 

contribution  of  the  bound  part  of  the  dielectric  function  COmes  in  a 

frequency  region  where  the  negative  contribution  of  the  free  electron  part 

of  the  dielectric  function  is  of  too  large  a  magnitude  for  to  reach 

zero.  It  is  the  aim  of  this  study  to  determine  the  optical  properties  of 

1  1 

copper  using  attenuated  total  internal  reflection  (ATR)  techniques  ’  and  to 
determine  if  a  surface  plasmon  resonance  exists  in  the  visible  region. 

Section  II  contains  the  experimental  details  of  sample  preparation  and 
surface  plasmon  resonance  measurements  and  analysis.  In  Section  III  ve  give 
the  results  of  the  optical  measurements,  namely  the  optical  properties  of 
copper  and  the  dispersion  of  its  surface  plasmon.  This  study  is  concluded 
vich  a  summary  in  Section  IV. 


-3- 


II.  EXPERIMENT 


A.  Sample  Preparation 

The  copper  films  were  prepared  by  chermal  evaporation  in  a  high  vacuum 
deposition  system.  All  films  were  prepared  on  the  same  day  that  the  experi¬ 
mental  measurements  were  performed,  so  as  to  reduce  the  effects  of  oxide  forming 
on  the  surface.  The  deposition  was  done  from  an  Al^O^  coated  evaporation  boat 
with  99.999"  pure  copper.  The  evaporation  was  performed  in  a  diffusion  pumped 
liquid  nitrogen  trapped  deposition  system  with  the  pressure  prior  to  growth  ap¬ 
proximately  2  x  10  ^  torr.  The  substrates  were  fused  silica  front  surface  mirror 
subscrates  (7.75  mm  diameter  *  4.00  am)  polished  to  better  than  1 / 10  wave 
flatness  at  0.5461  y  over  the  center  6  am  portion.  The  deposition  race  was 
150  Angstroms/ainute  and  was  monitored  during  evaporation  with  a  Sloan  Omni  11A 

rate  and  thickness  monitor.  The  total  thickness  of  the  Cu  films  was  3C00 
°  ’  Td 

Angstroms  as  measured  by  a  Sloan  Angstrometer "  .  The  prepared  films  had  a 
bright  metallic  luscer  and  visually  appeared  to  be  optically  smooth  and  homo¬ 
genous. 

3.  Experimental  Measurements 

t  ■? 

The  Otto  ATR  method*"  “  of  exciting  surface  plasnons  was  used  where  the 
surface  plasmon  was  for  an  air-Cu  interface.  The  ATR  measurements  were  per¬ 
formed  by  measuring  the  ratio  of  the  intensities  of  p-poiarized  light  (R_; 
to  thac  of  s-poiarized  light  (R^)  with  the  sample  pressed  against  the  base 
of  the  ATR  prism  and  normalized  to  the  reflectivity  with  the  sample  removed. 

The  prism  used  for  optical  measurements  was  a  (-5.0°  -  99.9°  . apex j  -  -o.l° 
light  flint  prism  which  was  flat  to  greater  than  a  tenth  of  3  wave  at  0.3-*bl  ... 
The  air  gap  was  obtained  by  using  either  evaporated  spacers  or  natural  iust 
The  uniformity  of  the  air  cap  was  adjusted  bv  observing  the 


particles . 
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uniformity  of  Che  surface  plasmon  resonance  of  copper  ac  angles  greater  than 
the  critical  angle  of  total  internal  reflection.  3eauti£ul  colors  of  red, 
green,  and  blue  could  be  visually  seen  for  p-polarized  light  depending  on  the 
viewing  angle  and  the  air  gap  thickness.  These  colors  could  be  seen  for  the 
largest  viewing  angle  possible,  therefore  indicating  the  wave-vector  divergence 
of  the  resonance  condition.  At  angles  less  than  the  critical  angle,  hewton 
interference  colors  were  not  visible  due  to  the  smallness  of  the  air  gap  used 
(d  a,  1000  A) .  The  reflection  spectra  were  taken  on  a  Perkin  Elmer  model  301 
far  infrared  spectrometer  modified  to  operate  in  the  single  beam  mode  in  the 
visible  with  the  addition  of  1440  lines/mm  grating,  appropriate  filters,  and 
use  of  tungsten  light  source  and  E.M.I.  photomultiplier  no.  9592B.  The  data 
collection  system  consisted  of  locking  in  on  the  13  cps  chopped  signal  with  a 
PAR  model  124A  Lock-In  Amplifier.  The  data  were  digitized  and  punched  on  paper 
tape  and  a  POP  11/20  computer  averaged  the  data  and  generated  a  plot  of  normal¬ 
ized  reflectivity  versus  wavelength.  A  digital  controller  was  used  to  sequence 
stepping  motor  changes  for  both  polarizations  and  for  stepping  the  wavelength 
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drive  through  the  visible  spectrum.  All  experimental  data  were  obtained  by 
fitting  Che  angle  of  incidence  of  the  input  lignc  and  by  varying  its  frequency; 
that  is,  the  complex  frequency  behavior  of  the  surface  plasma  resonance  was 
analyzed. 


C .  Optical  Constant  Determination 

The  ATR  experiments  were  performed  with  a  constant  air  gap  end  incident 

angle;  the  reflectivitv  ratio  (R  ,'R  )  was  measured  as  a  function  of  wave- 

p  s 

length.  All  optical  measurements  were  taken  at  room  temperature  and  are 

plotted  a3  normalized  reflectivity  ratio  lR  R  )  versus  wavelength,  as  shown 

P  ^ 

in  Figs.  1  and  2. 

To  theoretical  curve  fit  an  experimental  ATR  curve  the  wavelength 
dependent  optical  constants  n  and  k  need  cc  be  determined  along  with  the 

For  each  ATR  spectrum  the  wavelength 


wavelength  Independent  air  gap  1. 


dependent  normalized  reflectivity  ratio  (R  . R  )  was  measured  giving  us  ;r.e 
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set  of  values  for  two  unknown  wavelength  deoendent  values  n  and  k,'  and 
wavelength  independent  value  of  d.  Consequently,  it  is  impossible  to  determine 
the  optical  constants  and  the  air  gap  directly  from  the  exoerinentallv  obtained 
data.  The  method  that  was  chosen  was  to  use  the  experimencailv  determined 
optical  constant  k  obtained-irL^he  spectral  region  of  0.44  u  to  0.55  •_  bv 
L.  G.  Schultz. 4  tie  obtained  accurate  values  of  k  from  normal  incidence  trans¬ 
mission  measurements,  which  is  more  sensitive  for  determining  k  values  than 
reflectance  measurements.  The  spectral  region  of  0.44  u  to  0.55  u  was  chosen 
for  two  reasons,  the  first  being  that  the  relative  change  of  k  is  nearly  flat 
on  the  low  energy  side  and  slopes  off  slightly  on  the  high  energy  side;  the 
second  being  that  this  region  is  at  higher  energy  than  the  surface  plasnor. 
minimum  and  therefore  will  have  a  small  effect  on  the  optical  constants  oeter- 
mined  in  this  region. 

The  optical  constants  in  the  0.44  u  to  0.55  region  were  determined  by 

using  Schultz's"*  k  values  and  the  experimental!-.’  obtained  R  R  values  for  two 
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experimental  curves.  The  procedure  used  was  to  pick  an  arbitrary  air  vaiue  d, 
for  experimental  curve  1  and  then  to  calculate  from  the  R^  R  values  and  k 
values  a  sec  of  n  values.  Tith  this  set  of  n  and  k  values,  curve  1  was  com¬ 
puted  and  compared  to  the  exoerimentai  curve.  The  air  gap  for  the  second 
curve  d.,  was  varied  until  an  optimum  fit  was  obtained.  This  procedure  was  per¬ 
formed  several  times  by  varying  the  air  gap  i.  until  the  best  fit  for  curve  1 
was  obtained.  The  k  values  of  Schultz  were  accurate  tc  r  IT.  ana  therefore 
they  were  varied  within  this  range  to  optimize  the  fit.  Once  this  procedure 
was  completed,  a  set  of  n  and  k  values  was  determined  for  the  spectral  ranee 
0.44  u  to  0.55  u,  and  the  best  fit  3ir  gacs  i.  and  d,  were  determined.  The 
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final  scepoing  interval  for  d  was  r  5  Anestroms  whereas  the  n  values  were 
determined  to  0.0002  for  each  curve.  Once  the  air  gaps  a.  and  d,  were  determine 
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the  experimental  reflectance  data  could  be  used  directly  to  obtain  both  r.  and 
k  in  the  surface  plasraon  region.  The  procedure  for  this  region  vas  to  pick  a 
k  value  and  to  find  the  set  of  n  values  which  would  fit  the  experimental  curves. 
The  constraint  used  in  this  part  of  the  spectral  range  was  that  the  k  values 
should  vary/  smoothly  in  this  inter/al  and  that  the  error  be  displayed  in  the  n 
values.  Tne  reason  for  this  will  be  explained  in  tr.e  results  section.  Computer 
programs^  were  written  using  the  Fresnel  reflection  coefficients  and  the  three 
medium  multilayer  reflectance  formula  to  theoretically  computer  fit  the 
experimental  data.  The  computer  fits  are  shcvn  by  x's  on  the  same  graphs  as 
the  reflectivity  data  in  Figs.  1  and  2. 


III. 


RESULTS  OF  THE  OPTICAL  MEASUREMENTS  OF  THE  SURFACE  PLASMON  OF  COPPER 


A.  Optical  Properties 

The  AIR  determined  optical  constants  n  and  k  are  shown  in  Figs.  3  and  4, 
while  the  dielectric  functions  and  £2^'  are  shown  in  Figs.  5  and  6 

along  with  the  calculated  normal  incidence  ref. .activity  in  Fig.  7.  The  error 
in  determining  the  optical  constants  is  difficult  to  estimace.  The  average 
percentage  error  in  the  determined  index  of  refraction  as  shown  in  Fig.  3  was 
1.38%.  The  average  experimental  error  was  £  1%  which  gives  an  error  of 
-  3%  for  the  index  of  refraction  values,  while  an  error  of  2%  in  the  k 
values  yields  an  3.5%  error  in  n.  Combining  all  the  error  producing  terms 
yields  an  approximate  error  of  8%  for  the  n  values  and  4%  for  the  k  values. 

The  determined  optical  constants  are  quite  adequate  for  determining  the  opti¬ 
cal  properties  of  the  surface  olasmon  dispersion  which  is  the  main  purpose  of 
this  investigation.  3v  evaluating  the  imaginary  part  of  the  dielectric 
function  s?,  we  obtain  a  strong  absorption  edge  at  2.15  aV  which  is  attributed 
to  the  d-band  to  Fermi  surface  transitions  1.  -  L1  (E_)  and  L,(Q+)  - 
the  band  nomenclature  is  due  to  Mueller  and  Phillips.3  The  rapid  rise  in  t, 
is  due  to  the  flatness  of  the  d-band,  which  permits  a  large  number  of 
intarband  transitions  to  the  Fermi  surface  within  a  narrow  energy  range. 

The  first  peak  in  is  seen  at  2.48  eV  and  corresponds  to  the  transition 
L,(0-)  E_(L.,~)  .  These  values  agree  well  with  those  obtained  by  Nelkowsky 

J  T  *- 

and  3raunstein  which  were  2.18  eV  and  2.53  eV . 

3.  Dispersion  Analysis 

The  experimental  dispersion  of  the  surface  plasmon  resonance  can  be 
seen  by  noting  the  shift  of  the  minimum  for  two  different  incidence  angles 
in  Figs.  1  and  2.  The  ATR  experimental  minima  must  be  corrected  for  the 
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effects  or  the  prism.  These  effects  are  essentially  eliminated  by  using  only 
the  minima  obtained  when  che  minimum  reflectivity  for  o-poiarized  light  is 
greater  than  90%.  The  dispersion  curve  for  the  A7R  analysis  is  accomplished 
by  plotting  the  value  of  the  energy  at  the  minimum  of  each  AIR  curve  obtained 
for  different  incident  angles  (different  wavevectors  Sc..)  for  continuous  changes 

i 

in  the  frequency  of  the  incident  light.  The  AIR  computed^  dispersion  curve 
is  shown  in  Fig.  8  where  the  error  bars  signify  that  minima  were  determined 

a  a 

within  c  10  Angstroms  at  low  values  of  wavevector  and  :  1  Angstrom  at  high 

wavevector  points.  The  dispersion  curve  shows  the  characteristic  surface 

plasmon  resonance,  that  is,  ki  ;  *  •  at  'ji  «  a  .  The  computed  and  experimental 

ATR  spectra  had  less  than  10%  absorption  for  s-polarized  light  as  expected 

theroetic3lly .  The  resonance  value  for  the  plotted  wa-'e  vectors  is  at  2.163  eV . 

The  energy  wavevector  light  line  generated  bv  the  prism  for  the  experimental 

runs  of  Figs.  1  and  2  are  shown  on  che  dispersion  curve  plot  in  Fig.  8. 

The  experimental  ATR  minima  shown  in  Figs.  1  and  2  occur  at  lower  energy- 

due  to  che  prism  and  che  large  amount  of  overcoupling  C1-  75%  absorption  of 

incident  light) ;  che  minimum  is  partly  governed  by  the  resonance  condition 

C.  (uj)  *  -n  which  causes  the  resonance  to  shift  to  a  frecuer.cv  where  (c.)  is 
-  ?  ' 

more  negative.  The  location  of  the  surface  plasmon  resonance  is  shown  in  the 
plots  of  the  optical  constants  n  and  k  in  Figs.  3  and  1  and  on  the  dielectric 
function  plots  of  and  in  Figs.  5  and  6.  In  the  region  of  the  surface 
plasmon  resonance,  the  error  bars  for  determining  the  index  of  refraction  is 
large,  as  shown  in  Fig.  3.  This  error  is  mainly  due  to  the  effects  of  verv 
thin  oxide  layers  on  the  copper  surface  which  cause  the  resonance  condition 
to  shift.  This  is  the  reason  why  the  arbitrary  criteria  of  fixing  the  values 
of  k  to  lie  on  a  smooth  curve  was  chosen  as  explained  in  the  previous  section 
for  determining  the  optical  constants.  As  seen  from  Fig.  5,  the  surface  plasmon 
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resonance  occurs  ac  (*i)  =■  -■*.  7  which  is  a  large  deviation  from  the  norma* 

i. 

condition  e, (x)  =»  -1.  This  deviation  is  due  to  the  relatively  large  value  of 
£.,(:*))  =*  4.4  which  causes  the  resonance  to  be  damped  and  the  true  resonance 
to  occur  in  the  complex  frequency  plane  at  S(.~=ui  )  =•  -i.  It  should  also  be 

5 

noticed  that  the  resonance  occurs  on  the  3teep  absorption  edge  of  the  d-bar.d 
to  Fermi  level  transitions. 

Shown  in  Fig.  7  is  the  location  of  the  surface  plasmon  resonance  } 

s 

on  the  normal  incidence  reflectivity  curve.  This  resonance  should  be  barely 


visible  in  high  resolution  normal  incidence  reflectivity  spectra  on  copper 
which  has  a  rough  surface.  The  rough  surface  acts  as  a  grating  coupler  and 
permits  coupling  becveen  the  incident  transverse  light  and  longitudinal 
surface  plasmon  resonance. ' 

g 

Steel  has  done  a  study  on  the  characteristics  of  the  energy  loss  function 
(L)  for  Au  and  Ag.  He  used  a  Drude-lorencz  model  for  characterizing  the  di¬ 
electric  functions  of  Au  and  Ag  and  varied  the  free  electron  plasma  frequency 
and  relaxation  time  ana  studied  their  relationship  to  peaks  in  the  energy  loss 
function.  He  determined  that  "the  presence  of  the  absorption  band  causes  the 
structure  in  (L)  to  change  from  the  single  tree-electron  peak  to  twc  peaks: 
a  comparatively  weak  low-energy  peak  iLEF)  below  the  absorption  banc  and  also 
a  strong  high-energy  peak  (HE?)  above  X  " ;  where  is  the  bulk  plasma 
frequency.  The  surface  plasmon  resonance  of  copper  as  obtained  in  this  study 
is  equivalent  to  the  surface  plasmon  resonance  associated  with  the  LE?'  as 

O 

discussed  by  Steel"':  but  there  is  no  peak  in  the  energy  loss  function  at  these 
low  energies  where  the  energy  loss  function  for  the  surface  plasmon  resonance 
is  defined  at  1  >  -In[ 1/ Cs-1) ] .  The  , * -cable  reason  for  this  is  that  the 
damping  is  larger  in  Cu  than  it  is  for  both  Ag  and  Au  and  that  the  resonance 
is  truly  a  function  of  complex  frequency.  It  can  be  easilv  seen  that  the 


location  of  the  surface  plasmon  resonance  in  Fig.  7  is  in  the  region  of  the 
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steep  drop  in  normal  incidence  reflectivity.  This  reflectivity  drop  is 
normally  attributed  to  copper's  interband  transition  but  from  this  study  we 
can  see  that  reflectivity  drop  may  also  be  partially  caused  by  the  normal 
decrease  in  reflectivity  as  the  frequency  w  goes  through  the  surface  plasmon 
frequency  and  approaches  the  bulk  plasmon  frequency.  The  existence  of  the 
bulk  plasma  frequency  is  still  questionable  and,  if  it  exists,  it  will  occur 
in  the  complex  frequency  plane. 

It  was  of  interest  to  see  if  back-bending  of  the  complex  k  dispersion 
curve  occurred.  The  results  of  the  analysis  are  shown  in  Fig.  9  where  the 
real  part  of  the  complex  wavevector  (k^)  is  plocted  versus  energy.  The  plot 
of  the  complex— jj  behavior  is  also  shown.  It  can  be  seen  that  the  back-bending 
of  the  complex  k  dispersion  curve  occurs  when  the  complex-u  curve  diverges, 
therefore  having  the  desired  characteristics  of  the  surface  plasmon  resonance. 
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IV.  SUMMARY 

Experimental  results  show  chac  the  surface  plasmon  resonance  of  copper 
exists  in  the  visible  spectral  region  at  2.163  eV.  The  optical  properties 
of  copper  have  been  experimentally  determined  by  ATR  methods  and  used  to 
calculate  the  complex  frequency  dispersion  curve  which  shows  the  character¬ 
istic  wavevector  divergence  of  the  surface  plasmon  resonance,  while  the  complex- 
k  dispersion  curve  exhibits  the  characteristic  back-bending.  The  location 
of  the  resonance  is  on  the  steep  absorption  edges  of  £,  which  is  due  to  the 
d-band  to  Fermi  level  transitions.  These  transitions  cause  the  large  shift 
in  the  location  of  the  real  frequency  excitations  and  for  the  true  resonance 
to  occur  in  the  complex  frequency  plane.  The  color  of  the  noble  metal  of  copoer 
may  be  partially  attributed  to  the  steep  decrease  in  reflectivity  which  is 
caused  by  the  normal  decrease  in  reflectivity  as  the  excitation  frequency  m 
goes  through  the  surface  plasmon  frequency  and  approaches  the  bulk  plasmon 
frequency  and  to  copper's  numerous  interband  transitions.  The  existence  or 
the  bulk  plasma  frequency  is  still  questionable  and  if  it  exists  it  will  occur 
in  the  complex  frequency  plane. 
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.'..gure  -actions 


Fig.  1.  Attenuated  local  internal  reflection  spectrum  of  copper  plotted  as 
che  normalized  value  of  the  reflectivity  ratio  (R  'R  )  versus  wave- 

?  5 

length.  Solid  points  denote  the  experimental  data  with  a  solid  line 
drawn  through  them  for  clarity.  The  x's  are  che  compucer  fitted 
values.  Shown  in  the  lower  left  is  the  best  fit  air  gan  value  (d)  and 
the  angle  of  incidence  of  the  incoming  light  with  respect  to  the 
input  side  of  the  prism  as  shown  in  the  ATR  diagram  in  the  insert. 
Attenuated  total  internal  reflection  spectrum  of  copper  plotted  as 


•ig- 


Fig.  3. 
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Fig.  6. 


the  normalized  value  of  the  refleccivicv  ratio  (R  /R  )  versus  vave- 

?  s 

length.  Solid  points  denote  che  experimental  data  with  a  solid  line 
drawn  through  them  for  clarity.  The  x's  are  the  computer  fitted 
values.  Shown  in  the  lower  left  is  the  best  fit  air  gap  Value  i) 
and  the  angle  of  incidence  of  the  incoming  light  with  respect  to  the 
input  side  of  the  prism  as  shown  in  the  ..TR  diagram  in  the  insert. 

Plot  of  Che  real  part  of  che  refractive  index  of  copper  versus  wave¬ 
length  for  che  computer  fitted  data  of  Figs.  1  and  2.  Also  shown 
is  the  location  of  the  surface  piasmon  resonance  "•  )  . 

Plot  of  the  imaginary  part  of  the  retractive  index  of  copper  versus 
wavelength  for  the  computer  fitted  data  of  Figs.  1  ana  2.  Also 
shown  is  the  location  of  the  surface  piasmon  resonance  '■  ;  . 

Plot  of  the  real  part  of  the  dielectric  function  of  cooper  versus  wave¬ 
length  for  the  computer  fitted  data  of  Figs.  1  3nd  2.  Also  shown 
is  the  location  of  the  surface  piasmon  resonance 

s 

Plot  of  che  imaginary  part  of  the  dielectric  function  of  :cprer  versus 
wavelength  for  the  computer  fitted  data  of  Figs.  1  ana  2.  Vlso 


shown  is  the  location  of  the  surface  nlasmor.  resonance 
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Fig.  ?Io c  or  the  calculated  normal  Incidence  reflectivity  of  copper  using 

the  average  optical  constants  obtained  from  the  computer  fitting  of 

Figs,  i  and  2.  Also  shown  is  the  location  of  its  surface  plasmon 

resonance  (\  )  on  the  steep  reflectivity  edge. 
s 

Fig.  3.  Dispersion  curve  of  the  surface  plasmon  of  copper  (solid  line)  for 

vavevector  ki : .  F.ange  bars  indicate  that  the  minima  were  determined 

3  0 

within  z  10  A  for  low  vavevector  and  z  1  A  for  high  values  of  vave¬ 
vector.  Also  shown  in  the  energy-wavevector  light  line  generated  by 
the  prism  for  the  experimental  spectra  of  Figs.  1  and  2  where  the 
low  vavevector  line  corresponds  to  Fig.  1  and  the  high  vavevector 
line  to  Fig.  2.  The  angle  (AO)  is  with  respect  to  normal  to  the 
prism  base. 

Fig.  9.  Complex-'jj  and  complex-*  dispersion  curves  of  the  surface  plasmon  of 
copper  along  with  3  plot  of  the  energv-wavevector  relation  for  light 


in  air. 
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